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Abstract 

This paper aims to apply artificial neural network 
(ANN) based phase shifting transformer (PST) to 
enhance multi-machine power system transient 
stability. A new model of the phase shifting 
transformer based artificial neural network is 
proposed. The proposed model depends mainly on 
the function of the phase shifting transformer. The 
output machine voltages deviations, speed 
deviations and angles deviations are used as an 
input signal to the ANN model to provide phase 
advance for each machine voltage of the studied 
three machine six bus interconnected power 
system. The ANN model consists of input layer, 
hidden layer and output layer. The input layer 
consists of 6 input signals. The output layer 
consists of one signal which is the phase advance. 
The ANN offline training is made first to initialize 
the weights and bias matrix values. Hence, adaptive 
function of ANN is used as online ANN application 
to the studied multimachine power system to 
modify the weights and bias matrix according to 
the system dynamic performance. Different fault 
locations were considered to judge the 
effectiveness of the proposed ANN based PST. The 
time simulation results prove that the proposed 
ANN model based PST is very effective in 
improving the power system transient stability in 
case of severe disturbance such as unrepeated 
three-phase short circuit fault. A Comparative 
study between the conventional PST and ANN 
based PST proves the superiority of the proposed 
model. The studied power system is modeled and 
solved using the MATLAB software package. 

Keywords: Transient Stability enhancement - 
Multi-machine Power System - Phase shifting 
transformer - Artificial neural network model. 



l.Introduction 

Specially designed power system devices can 
control the flow of the active and/or reactive power 
in interconnected power systems alleviating line 
overloads and changing parallel line load sharing. 
Phase shifting transformer (PST) was the only 
device capable of controlling both magnitude and 
direction of the active power flow. Nowadays, the 
traditional PST is still the most widespread 
technology!!]. Phase shifting transformer (PST) can 
provide power flow control in power systems, by 
varying the phase angle between two buses. The 
variation of phase angle is carried out by adding a 
regulated quadrature voltage to the source line to 
neutral voltage. The use of PST is better than 
generation rescheduling from the economic point 
of view[2]. The PST can also be used for power 
system stability enhancement, and to increase the 
transmission capacity of the electric power 
network[3-9]. PST is considered the principle 
control device that can direct power flow in a 
specific part of the electric network[10]. The Phase 
shifting transformer is an electric device from the 
FACTs (Flexible AC Transmission system) devices[9]. 
A Proper method to determine the optimal 
allocation pattern of multiple phase shifting 
transformers for unwinding parallel flows has been 
investigated and studied[ll]. The specification, 
design and protection of phase shifting 

transformers for enhanced interconnected real 
power system have been studied[12]. Due to 
growth in electric power networks to a large-scale 
complicated system, an increase in control system 
is imposed to meet the global stability 

improvement. Artificial neural network can be 
tuned to change system conditions and to provide 
better control signal than a conventional control 
approach [13]. The applications of artificial neural 
networks (ANN) in control, and modeling to 
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electric power systems have been investigated for 
different electric devices and different operating 
power systems parameters[14-20]. The ANN 
control approach is widely applied to all flexible AC 
transmission system devices (FACTS). Also, The 
ANN is used to model the nonlinear loads, wind 
turbine, and microturbine controller[21-23]. 

This paper proposes an application of Artificial 
Neural Networks (ANN) based phase shifting 
transformer (PST) to improve the transient stability 
of interconnected multimachine power systems. 
The PST is modeled based ANN approach. The 
proposed ANN model can provide the function of 
the PST and be applied successfully online with the 
studied multimachine power system. The proposed 
model of the PST based ANN is simulated with the 
studied 3-machine 6-bus power system to judge 
the effect of inserting the proposed ANN based 
PST in improving the transient stability in case of 
severe disturbance. The ANN is used first for offline 
training to obtain the initial values of the weights 
and bias matrix. Hence adaptive ANN function is 
applied with the studied system, which provides 
online ANN based PST. Different fault locations 
were considered in this study to judge the effect of 
the proposed approach. The time simulation 
results proves the effectiveness of insertion of PST 
based ANN in enhancing the transient stability of 
the power system in case of severe disturbance 
such as unrepeated three phase short circuit fault. 
The comparative study proves the superiority of 
the proposed ANN based PST over the 
conventional PST. The study considers that the 
voltage is constant and just change the phase 
angle of the phase shift transformer. 



2. Studied Power System Model 

The multimachine power system used to evaluate 
the effectiveness of installing Phase Shift 
Transformer based ANN in enhancing the transient 
stability is shown in figure 1[24]. 



A reduced order model is used for this study. An 
energy function of an electric power system is 
constructed as a first integral of motion(swing 
equation). Using the system transient state, the 
following swing equation can be written as 
follows[25]: 



M dA® _ p p 

dt m e dt 



( 1 ) 



dS_ 

dt 



= A co 



( 2 ) 



Where the generator is modeled in the standard 
simplified way: A constant mechanical power, an 
electromotive force of constant modulus behind a 
transient reactance, and lossless transmission 
system. The modification of a real electric power 
transfer is given as: 




power 



P e = -^-^-sin(<^ +a) (3) 

X 

3. Phase Shifting Transformer 

The phase shift transformer controls the active 
power flow in interconnected power systems by 
varying the phase angle between two buses. The 
variation of this phase angle is carried out by 
adding a regulated quadrature voltage to the 
source line-to-neutral voltage. PSTs have different 
forms of characteristics, which are: 

• Direct PST's are based on one 3-phase core. The 
phase shift is obtained by connecting the winding 
in a suitable manner. 
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• Indirect PST's are based on a construction with 
two separate transformers. One variable tap exciter 
to regulate the amplitude of the quadrature 
voltage, and series transformer to inject the 
quadrature voltage in the right phase. 

• Unsymmetrical PST's create an output voltage 
with an altered phase angle and amplitude 
compared to the input voltage. 

• Symmetrical PST's can create an output voltage 
with an altered phase angle compared to the input 
voltage, but with the same amplitude. 

The combination of the previous characteristics can 
results in 4 categories of PST's. In this study there is 
only one used and will be discussed in the 
following paragraph. 




3.1 Direct Symmetrical PST , s[26] 

By additional tap changer, the direct 
unsymmetrical PST can be made symmetrical, the 
schematic diagram is shown in figure 2. The 
advantages are increase in maximum attainable 
angles and that the amplitude voltage remains 
unchanged. The relation between the angle a and 
the quadrature voltage is still nonlinear and can be 
derived from the phasor diagram shown in figure 



3. 




direct, 



3.2. Phase Shifting Transformer Model 

The phase shift transformer connected between p 
and q buses can be represented by the equivalent 
circuit shown in figure 4. [27]. It consists of an 
admittance in series with an ideal transformer 

having a complex turns ratio n = nZa. The phase 
angle step size may not be equal at different tap 
positions. However, equal step size is usually used 
in power flow and transient stability programs. 




Figure 4 Phase shifting transformer 
representation 



The terminal voltages are related by the following 
equation: 



— — = nZa = n ( cos a + j sin a) - o s + j b s (4) 



Where a is the phase shift from bus p to bus q, it is 

positive when v p leads v q . Since there is no 
power loss in ideal transformer, 



. * 

v p 1 p 



. * 

' V q l s 



Therefore, the transformer current at bus p is: 

: i:y 



l p — 



h = - 



- A a s - A 



-{Vq V.y ) 



(5) 



( 6 ) 
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Substitute for v q from the previous equation we 
can get: 

Y 

ip = 2 6 2 t V P~( a s + J b s) V s] (7) 

a s +b S 

Substitute for i p the value of i s can be obtained 
as the following equation: 

is — j[~v P +(a s + jb s )v s ] (8) 

a s + b s 

From the previous equation, the phase shifter 
terminal voltages and currents can be arranged in 
the following matrix equation as: 







1 

1 


j- 


1 P 




a2 s + b s a s ~ J b s 


V P 






Y. 


*v 


h. 




Y e 

L a s + J b s 


_V,_ 



It is clear that the previous admittance matrix is 
unsymmetrical, so that the transfer admittance 
from bus p to s is not equal the transfer admittance 
from s to p. The previous equation incorporating 
with the power system model and solved using 
Matlab software. 



3. Output signals from the hidden nodes result 
from the input signals passing through their 
tansigmoid activation function. 

4. Hidden layer output signals are sent to output 
nodes through weighted connections between the 
hidden nodes and the output nodes. 

5. The ANN output signal is obtained using 
tansigmoid transfer function of the output node. 




Figure 5 Neural Network Model of PST 



4. Neural Network Model of Phase 
shifting Transformer 

To demonstrate the neural network model of the 
phase shifting transformer, different operating 
conditions, different fault duration and different 
fault location are considered to obtain the studied 
power system performance data. The data 
obtained are machine terminal voltages 
magnitude, angle deviations and speed deviation. 
Hence the neural network construction consists of 
input, hidden and output layers. The input signals 
are voltages, speed deviation, angle deviations and 
the output signal is the phase shift angle. 

4.1 Initialize the Neural Network Data 
The proposed ANN based PST consists of three 
layers, named input layers, hidden layer and output 
layer. The input layer has six nodes. After try and 
error the best number of hidden layer found as 2 
nodes with nonlinear tansigmoid function. The 
output layer consists of one node. The ANN 
structure is shown in figure 5. The operation of 
such ANN can be described as follows[28]. 

1. The input layer nodes receive the data from the 
studied power system. 

2. The output signals from the input nodes pass to 
hidden nodes via the weighed links which can be 
gained or attenuate. 



The ANN desired output signal is the phase shift 
angle, where the voltage is considered to be 
constant value of 1.00 p.u.. The statistical data of 
the ANN training are as follows: 

No. of training data = 2000 
No. of iterations = 8000 
Sum Squared Error = 0.0001 
Learning rate = 0.001 
Momentum constant = 0.96 

The resulting weight matrix between the hidden 
and input neurons is given by: 

[0.9087 0.8024 -0.6704 0.3963 -0.9685 -0.9532' 

W, = 

1 [0.7601 -0.9599 0.5067 0.3897 -0.6077 0.6483 

The hidden nodes biases are: 

[- 0 . 8688 ] 

B x = 

1 -0.3811 



The connection weights between the output 
neurons and the hidden neurons are: 

W 2 = [-0.5398 0.1837] 

The bias of the output neurons is : 

B 2 = [ -0.3144] 

4.2 Online Neural Network based PST 



The input layers receives Acj(k), Auj(k-l), AS(k), 
A5 t (k-1), Av t (k), Avt(k-l). The back propagation 
learning algorithm is an interactive method 
employing the gradient decent algorithm for 
minimizing the mean square error between the 
actual output and the target for each pattern in the 
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training set. The back propagation algorithm can 
be discussed in the following steps: 

1. Initializing the network bias and weight matrix 
by any values or zeros. 

2. Calculating the actual output of the network 
layers. 

3. Updating the network weights in a recursive 
algorithm starting from the output layer and 
working backward to the first hidden layer. The 



weights and biases are updated as: 

Wy( t+1) = Wy{t) + A Wy (t) (10) 

0/t+l) = 0 j(X) + A 0 j(t) (11) 

Where: 

A wg (t) = p5j P O jp + pA Wy{ t-1 (12) 

A 0, (t) = p5j P (t) + PA0/t-l) (13) 

4. If j is an output layer, then the error will be: 

5 jP = F^j* (Tjp-Ojp) (14) 

5. If j is a hidden layer then the error will be : 

5 JP = F (a jp )* I S jp w LJ (15) 

Where: 



A w-y (t) The weight change from a neuron in the i th 
layer to the neuron in the j th layer 
Ojp The network output in the j th unit in the output 
layer for a pattern p. 

Tj P Target value in the training pattern 

|i The learning rate 

w-y The weights from layer i to layer j 

P The momentum factor 

t The iteration number 

Qj The bias for the layer j 

F The sigmoid activation function. 

6. Check the error goal. If the condition is valid, 
keep the recent data of the control system, 
otherwise start from the beginning. 

The previous process is repeated independently for 
each PST connected directly to each machine. 

5. Results and Discussions: 

The System under study is shown in figure 1. The 
parameters of the system are given in Ref. [24]. 
Firstly, the load flow study is made to obtain the 
bus voltages and angles at different operating 
conditions for the transient stability study. The 
conventional PST is directly connected to each 
machine, immediately after the three-phase short 
circuit fault is considered for different intervals and 
the system is solved to obtain the system 
performance. The data are saved to be used for the 
offline training of the proposed ANN. The weights 
and bias matrix obtained from the offline training 
are considered as an initial value of the proposed 
adaptive neural network. The ANN is updated 
online with the studied system to stabilize the 



system under study. The flowchart describing the 
solution of the system understudy is shown in 
figure 6. 




Figure 6 Flowchart of the time simulation of 
the studied power system with ANN based 



To judge the ability of the proposed artificial neural 
network based phase shifting transformer to 
improve the transient stability of the studied 3- 
machine six bus power system in case of severe 
disturbance. A Comparative study between the 
dynamic performance of the studied multimachine 
power system equipped by conventional PST and 
ANN based PST is obtained. Also the dynamic 
performance of the studied system without 
installing PST is given. 
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Fig.7 Studied system performance due 3-ph. s.c. 
fault at bus # 1 without installing PST 



5.1. System performance without 
installing PST 

Figures 7, 9 and 11 shows the studied system 
performance when 3-phase to ground short circuit 
fault is applied to bus No. 1, 3 and 6 without 
installing phase shifting transformer. It is very clear 
that the system is unstable. 




Fig. 8 Studied system performance due 3-ph. s.c. 
fault at bus # 1 with installing PST and ANN based 
PST 




Fig. 9 Studied system performance due 3-ph. s.c. 
fault at bus # 3 without installing PST 

5.2. System performance with installing 
PST and ANN based PST 

Figures 7, 9 and 11 shows the studied system 
performance when 3-phase to ground short circuit 
fault is applied to bus No. 1, 3 and 6 with installing 
conventional phase shifting transformer, and 
artificial neural network based phase shift 
transformer. It is very clear that the system is stable 
with both conventional PST and ANN PST. The 
comparative study shows that the ANN based PST 
has better performance than the conventional PST 
in damping the system oscillations very fast, with 
less overshoot and undershoot. Also the settling 
time with installing ANN based PST is very short 
compared with conventional PST. 




Fig. 10 Studied system performance due 3-ph. s.c. 
fault at bus # 3 with installing PST and ANN based 
PST 
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Fig. 11 Studied system performance due 3-ph. s.c. 
fault at bus # 6 without installing PST 




f\ 




1 1 — 


PST 

. — ANN based PST 

1 1 1 1 1 1 1 



) ' ' ' ' ' ' ' ' ' 

□ 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Time in sec. 

Figure 12 Studied system performance due 3-ph. 
s.c. fault at bus # 6 with installing PST and ANN 
based PST 

6. Conclusions 

This study proposes a novel artificial neural 
network model for the phase shifting transformer. 
The studied system performance under different 
operating conditions, different fault locations and 
different fault duration time is considered to 
provide enough data for the neural network 
training, to ensure a good result of the network 
parameters. The proposed ANN is constructed and 
the data is trained to obtain the network 
parameters. Hence adaptive neural network 
function is applied. The time simulation proved 
that the proposed ANN based PST has better 
performance than that conventional PST with 
better control quality, less overshoot and under 
shoot. 
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Abstract 

The aim of credit scoring system is to model or 
predict the probability that consumers with certain 
characteristics are to be considered as potential risks. 
Until recently, the decision to grant credit was based 
on human judgment to assess the risk of failure to 
pay. From the literature it has been found that data 
mining models such as support vector machine 
(SVM), multilayer perceptron neural network 
(MLPNN) and decision tree (DT) can help decision 
makers to improve in this domain. SVM, MLPNN and 
DT with their remarkable ability to derive meaning 
from complicated or imprecise data, can be used to 
extract patterns and detect trends that are too complex 
to be noticed by either humans or other conventional 
techniques. This paper evaluates the performance of 
SVM with polynomial kernel (POLY-SVM), MLPNN 
with pruning parameters and DT with Chi-squared 
automatic interaction detection (CHAID) algorithm in 
the prediction of consumer’s credibility. The dataset 
contains information related to the sociodemographic 
and financial characteristics of consumers from real 
world credit bank. Known sets of risky and 
creditworthy consumer data were used to train the 
models to categorize new cases. The purpose was to 
determine an optimum classification model with high 
predicting accuracy for this problem. Experimental 
results demonstrate the effectiveness of CHAID-DT 
model in predicting credit scoring with higher 
accuracies than other techniques. 

Keywords: Credit scoring, data mining, support 
vector machine, neural network, decision tree. 

1. Introduction 

Normally, consumers who applied for a loan or a 
credit, have to fill out a questionnaire that asks for 
relevant financial and personal information. This 
information is used by a bank (or a financial 
institution) as the basis for its decision as to whether 
to lend money or not. Until recently, the decision to 
grant credit was based on human judgment to assess 
the risk of failure to pay [3]. Actually, lending to 
consumers exploded over the last two decades. 



Consumer credit has had one of the highest growth 
rates in any sector of the business. Consumers use 
credit to acquire goods and services now and then pay 
for them later. Many credit applicants are, in fact, 
good credit risks, but some are not. The risk for 
financial institutions comes from not knowing how to 
distinguish the creditworthy applicants from the bad 
ones. There is a need to use more formal and objective 
methods (generally known as credit scoring) to help 
credit providers decide whether to grant credit to an 
applicant. Credit scoring is the set of decision models 
and their underlying methods that help financial 
institutions in the granting of consumer credit [4]. The 
goal of these methods is to predict credit applicants 
into two classes: the “good credit” or creditworthy 
class that is liable to reimburse the financial 
obligation and the “bad credit” or risky class that 
should be denied credit due to the high probability of 
defaulting on the financial obligation. The prediction 
is dependent on sociodemographic characteristics of 
the consumer (such as age, education level, 
occupation and income), the repayment performance 
on previous loans and the type of loan. From the 
literature it has been found that data mining models 
such as SVM, MLPNN and DT can help decision 
maker to improve in this domain [5-7]. They have 
remarkable ability to derive meaning from 
complicated or imprecise data, can be used to extract 
patterns and able to detect trends that are too complex 
to be noticed by either humans or other conventional 
techniques. This paper evaluates the performances of 
SVM with polynomial kernel, MLPNN with pruning 
parameters and DT with Chi- squared automatic 
interaction detection in the prediction of consumer’s 
credibility. The three models have been trained on 
data samples of real world credit dataset from the 
University of California Irvine (UCI) Machine 
Learning Repository [1]. The remainder of the paper 
is organized as follows: Section (2) focuses on the 
sociodemographic and financial characteristics of 
consumers and reviews the recent published papers in 
the credit scoring subject. Section (3) reviews the 
three classification models and presents their 
parameters that need to be optimized. Section 
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demonstrates the experimental results to show the 
effectiveness of each model. Finally, conclusions are 
introduced in Section (5). 

2. Background 

2.1. Credit card dataset 

The Australian credit dataset used in this dataset is 
obtained from the UCI Machine Learning Repository 
[1]. The dataset consists of I I I I samples □ 3D7 
correspond to creditworthy applicants and 3 EB 



correspond to risky ones to whom credit requests 
should be refused. Each sample contains 1 5 attributes 
of which If attributes are continuous (range) and the 
remaining attributes are categorical (nominal). In 
order to preserve the confidentiality of the data, the 
names and values of the attributes were replaced by 
meaningless identifiers. Table 1 shows the attributes 
of the dataset with the number of missing values for 
each one. The missing values are replaces with the 
mode or the mean of each attribute. 



Table 1. Contents of the Australian credit dataset 



Attribute 


Type 


Values 


Number of 
missing values 


□Mode Mean for 
attributes with 
missing values 


Al 


Flag 


b, a. 


12 


b 


A2 


Range 


Continuous. 


12 


31.57 


A3 


Range 


Continuous. 






A4 


Nominal 


u, y, l,t 


□ 


u 


A5 


Nominal 


g> P,gg 


□ 


g 


AD 


Nominal 


c, d, cc, i, j, k, m, r, q, w, x, e, aa, ff. 


D 


c 


A7 


Nominal 


v, h, bb, j, n, z, dd, ff, o. 


it 


V 


AD 


Range 


Continuous. 






AD 


Flag 


t, f. 






Al | 


Flag 


t, f. 






All 


Range 


Continuous. 






A12 


Flag 


t, f. 






A13 


Nominal 


g, P, s. 






A14 


Range 


Continuous. 


13 


1 :4. 2 


A15 


Range 


Continuous. 






Class 


Flag 









Qnissing values are replaced with the mode or the mean of the attribute. 



2.2. Previous works 

The logistic regression has been for long time the 
most used statistical model in assessing the credit risk 
[4, Cj* Afterward, a number of quantitative models are 
being used also including linear discriminant analysis 
(LDA) [□] artificial neural network (ANN) [ID], 
ensemble of neural networks [11], decision tree [12], 
genetic programming [13] and support vector machine 
[14]. Two recent reviews of these models published in 
2 I 7 and 2D1Q can be found in [Ej and [15] 
respectively. Actually, the issue of credit scoring is 
attracting several scientists □ search for “credit and 
scoring” on sciencedirect library results in 14 □ studies 
in anuary 21 and 72 studies published in 2EDQ In 
the Egyptian literature, there are two remarkable 
papers which have been recently published in 2l I I I 
and 2 DID by Hussein Abdo et al. In [ID] they 
concluded that neural network models gave a better 
average correct classification rate than conventional 
techniques such as, discriminant analysis and logistic 
regression. In [13], they investigated the ability of 
genetic programming in the analysis of credit scoring 
models in the Egyptian public sector banks. 



3. Classification models 
3.1. Decision tree 

DT models are powerful classification algorithms. 
They are becoming increasingly more popular with 
the growth of data mining applications [5, □]. As the 
name implies, this model recursively splits data 
samples into branches to construct a tree structure for 
the purpose of improving the prediction accuracy. 
Each tree node is either a leaf node or decision node. 
All decision nodes have splits, testing the values of 
some functions of data attributes. Each branch from 
the decision node corresponds to a different outcome 
of the test. Each leaf node has a class label attached to 
it. Figure 1 shows an illustrating example for binary 
tree where each decision node has two splits only. 
CHAID-DT model allows multiple splits on a 
predictive attribute. This model relies on the Chi- 
2 

square X test to determine the best split at each step. 
CHAID-DT only accepts nominal or ordinal 
categorical predictive attributes [ID]. If hen predictors 
are continuous, they are transformed into ordinal ones. 
Ordinal attribute is ordered set with intrinsic ranking. 
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Root Node 




Figure 1. Example of a binary decision tree. 



The CHAID-DT modeling algorithm is as follows [ □] : 

1. Dinning the continuous attribute (if exists) to 
create a set of categories, where each category is 
a subrange along the entire range of the attribute. 
This binning operation permits CHAID-DT 
model to accept both categorical and continuous 
inputs, although it internally works only with 
categorical ones. 

2. Analyzing the categories of each attribute to 
determine which ones can be merged safely to 
reduce the number of categories. 

3. Computing the adjusted p value for the merged 
categories by applying Donferroni adjustments. 

4. Searching for the split point with the smallest 
adjusted p value (probability value, which can 
be related to significance) to find the best split. 

In step 2, the algorithm merges values that are judged 
to be statistically homogeneous (similar) with respect 
to the target attribute and maintains all other values 
that are heterogeneous (dissimilar). If the p value is 
greater than specified parameter a mer g then the 
algorithm merges the pair of categories into a single 
one. The value of a mer g must be greater than Dand less 
than or equal to 1. To prevent any merging of 
categories, specify a value of 1. In step 4, each 
predictive attribute is evaluated for its association 
with the target attribute, based on the adjusted p value 
of the statistical test of association. The predictive 
attribute with the strongest association, indicated by 
the smallest adjusted p value , is compared to a pre- 
specified split threshold a sp i it . If the adjusted p -\ alue is 
less than or equal to a sp i it , that attribute is selected as 
the split attribute for the current node. After the split 
is applied to the current node, the child nodes are 
examined to see if they warrant splitting by applying 
the merge Split process to each in turn. Processing 
proceeds recursively until one or more stopping rules 



are triggered for every unsplit node, and no further 
splits can be made. In this study, the target attribute is 
of categorical type (creditworthy or risky). The 
Likelihood ratio has been used to compute the chi- 
square statistic. The algorithm forms a contingency 
(count) table using the classes of the target attribute 
as columns and the categories of the predictive 
attribute as rows. The expected cell frequencies 
under the null hypothesis of independence are 
estimated. The observed cell frequencies and the 
expected cell ones are used to calculate the chi- 
squared statistic and the p value. 



=ZZ ln 



f \ 

n, 



=i i = i 



m 

V 1 J 



( 1 ) 



where Tl 



i ^ n ( n - ^ n ) 



is the 



observed cell frequency and m. is the expected cell 



frequency for cell ( n = /, n — ) , and the p 
value is computed as follows: 

P = Prfc> 2 ) ( 2 ) 

The CHAID-DT model is fast, builds “wider” 
decision trees as it is not constrained to make binary 
splits Quaking it very popular in market research. The 
model can be conveniently summarized in a simple 
two-way contingency table, with multiple categories 
for each variable. However, this algorithm requires 
larger quantities of data to get dependable results. 



3.2. Neural network 

Multilayer perceptron neural network (MLPNN) 
with back-propagation is the most popular 
artificial neural network architecture [17,1 □]. 
The MLPNN is known to be a powerful function 
approximator for prediction and classification 
problems. It is debatably the most commonly 
used and well-studied ANN architecture. Figure 
2 shows structure of MLPNN which is organized 
into layers of neurons input, output and hidden 
layers. There is at least one hidden layer, where 
the actual computations of the network are 
processed. Each neuron in the hidden layer sums 
its input attributes t after multiplying them by 
the strengths of the respective connection 
weights Wjj and computes its output using 
activation function (AF) of this sum. AF may 
range from a simple threshold function, or a 
sigmoidal, hyperbolic tangent, or radial basis 
function. 



i Wi i (3) 

where is the activation function. 

□ack-propagation (DP) is a common training 
technique for MLPNN. The available dataset is 
normally divided into training and test subsets. 
□P works by presenting each input sample to the 
network where the estimated output is computed 
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by performing weighted sums and transfer 
functions. The sum of squared differences 
between the desired and actual values of the 
output neurons is defined as 

= |l( „ - ) 2 (4) 

where d is the desired value of output neuron and 
is the actual output of that neuron. 

□ eights w t in Equation (3), are adjusted to 
reduce the error of Equation (4) as fast, 

quickly as possible. DP applies a weight 
correction to reduce the difference between the 
network outputs and the desired onesD i.e., the 
neural network can learn, and can thus reduce 
the future errors. The performance of MLPNN 
depends on network parameters, the network 
weights and the type of transfer functions used 

[iq. 



Creditworthy 



Risky 



J 

Input Hidden Output 

Layer Layer Layer 

Figure 2. The structure of multilayer perceptron 
neural network. 

□ hen using an MLPNN, three important issues need 
to be addressed □ the selection of data samples for 
network training, the selection of an appropriate and 
efficient training algorithm and determination of 
network size. New algorithms for data portioning [2q 
and efficient training with faster convergence 
properties and less computational requirements are 
being developed [21]. However, the third issue is a 
more difficult problem to solve. It is necessary to find 
a network structure small enough to meet certain 
performance specifications. Pruning methods for 
improving the input- side redundant connections were 
also developed that resulted in smaller networks 
without degrading or compromising their performance 
[ 22 ]. 

3.3. Support vector machine 

The support vector machine is a supervised machine 
learning technique, which is based on the statistical 
learning theory. It was firstly proposed by Cortes and 
Vapnik from his original work on structural risk 
minimization in [24] and modified by Vapnik in [25]. 
The algorithm of SVM is able to create a complex 




decision boundary between two classes with good 
classification ability. The basic idea is to map the data 
into a higher dimensional feature space via a 
nonlinear mapping kernel function chosen a priori 
(Figure 3), and constructs a hyperplane, which splits 
class members from non-members (Figure 4). 





Input space 



Feature space 



Figure 3. Mapping kernel. 



upport 



Optimal 

Hyperplane 



Figure 4. Optimal hyperplane. 

SVMs introduce the concept of Tnargin’ on either 
side of a hyperplane that separates the two classes. 
Maximizing the margins and thus creating the largest 
possible distance between the separating hyperplane 
and the samples on either side, is proven to reduce an 
upper bound on the expected generalization error. 
SVM may be considered as a linear classifier in the 
feature space. On the other side it becomes a 
nonlinear classifier as a result of the nonlinear 
mapping from the input space to the feature one [2[tf 
27]. For linearly separable classes, SVM divides these 
classes by finding the optimal (with maximum 
margin) separating hyperplane. Optimal hyperplane 
can be found by solving a convex quadratic 
programming (IOP) problem [5]. Once the optimum 
separating hyperplane is found, data samples that lie 
on its margins are known as support vectors. The 
solution to this optimization problem is a global one. 
For linearly decision space, suppose the training 
subset consists of n samples ( b i)*D ,( n ), E p 

and EZf]l,-l i.e. the data contains only two classes. 
The separating hyperplane can be written as 

(,) = w« + (5) 

where the vector w and constant are learned from 

a training subset of linearly separable samples. 

The solution of SVM is equivalent to solve a linearly 
constrained quadratic programming problem as 
Equation (Q for both targets §1 1 and □!. 
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i = w t + >1, idRn ( ) 

As mentioned before, samples that provide the above 
formula in case of equality are referred as support 
vectors. SVM classifies any new sample using these 
support vectors. On the other hand, margins of the 
hyperplane follow the subsequent inequality 

jX . , ( '^l', iU[,a,n. ( 7 ) 

The norm of the w has to be minimized in order to 
maximize the margin F . In order to lessen the number 
of solutions to the norm of w, the following equation 
is assumed 

rx H|=i (□) 

Then the algorithm tries to minimize the value of 

lC2||w|| subject to the condition (E). In the case of 

non-separable samples, slack parameters £ are added 
into the condition as follows. 

,.(W ,.+ Vi (q{ 

And the value that we want to minimize becomes: 

c^+iajwf. (iq 

i = 1 

C is a regularization parameter to determine the level 
of tolerance the model has, with larger C values 
allowing larger deviations from the optimal solution. 
This parameter is optimized to balance the 
classification error with the complexity of the model. 
There is a family of kernel functions that may be used 
to map input space into feature space. They range 
from simple linear and polynomial mappings to 
sigmoids and radial basis functions. Once a 
hyperplane has been created, the kernel function is 
used to map new samples into the feature space for 
classification. This mapping technique makes SVM 
dimensionally independent, whereas other machine 
learning techniques are not. In this study, we used the 
polynomial function kernel to map the input space 
into the higher dimensional feature space. Polynomial 
kernels can be controlled by adjusting the complexity 
of the mapping d , the coefficient r and weight value 
y in the kernel function (Equation 11) [2D]. 

( ) = h« • +r)\ym (11) 

where y , r and d are the kernel parameters 
Consequently, there are four parameters need to be 
optimized in the polynomial kernel-SVM model Q 
C, y , r and d. 

4. Experimental results 

Figure 5 shows components of the proposed stream 
for credit scoring. The stream is implemented in SPSS 
Clementine data mining workbench using Intel core 2 
Dup CPU with 2.ED DHz. Clementine uses 
client Server architecture to distribute requests for 
resource-intensive operations to powerful server 
software, resulting in faster performance on larger 



datasets [2]. The software offers many modeling 
techniques, such as prediction, classification, 
segmentation, and association detection algorithms. 
The Australian credit dataset contains :eD 3 of 
sociodemographic and financial data records each 
with 15 attributes and 2 target classes; creditworthy 
and risky. The whole dataset is divided for training 
the models, and test them. The training set is used to 
estimate the model parameters, while the test one is 
used to independently assess these models. These 
models are applied again to the entire dataset and to 
any new data. 

Credit cards dataset node is connected directly to an 
Excel file that contains the source data. The dataset 
was explored for incorrect, inconsistent or missing 
data. The missing values are replaced with a mode for 
nominal attribute or a mean for range one. The range 
attributes are normalized before the analysis using 
MLP and SVM models. On the other side, they are 
kept unchanged for the analysis of DT model. 

Data audit node provides a comprehensive first look 
at the attribute values in an easy-to-read matrix that 
can be sorted and used to generate full-size graphs. 
The audit node computes different descriptive 
statistics and histogram of each attribute in the 
dataset. 

Type node specifies the field metadata and properties 
that are important for modeling and other works in 
Clementine. These properties include specifying a 
usage type, setting options for handling missing 
values, as well as setting the role of an attribute for 
modeling purposes □ input or output. As previously 
stated, the first 15 attributes in Table 1 are defined as 
input (predictive) attributes and the credibility of a 
consumer is defined as target class. Predictive 
attributes are of mixed types □nominal, flag or range, 
while the target class is of type flag. 

Partition node is used to generate a partition field that 
splits the dataset into separate subsets for the training 
and test the models. In this study the dataset was 
partitioned by the ratio 7Q3CO for training and test 
subsets respectively. 

CHAID classifier node is trained using the 
significance of a statistical test as a criterion. The 
modeling algorithm evaluates all of the values of each 
potential predictive attribute. The chi-square statistic 
is computed using the likelihood ratio to merge values 
that are judged to be statistically homogeneous with 
respect to the target attribute. After merging similar 
categories, the algorithm selects the best predictive 
attribute to form the first branch in the decision tree, 
such that each child node is made of a group of 
homogeneous values of the selected attribute. The 
best values of (Verge and a sp i it are found to be □ 1 5 for 
both and the maximum allowable depth is set to 5. 
The resulting accuracies for training and testing are 
rTi nxi and J|T7CD respectively. CHAID-DT model 
is very fastDit takes below one second to build the 
model and easy to interpret Dit can be expressed as a 
set of rules (Figure Q. 
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Data Audit 



OH AID DT MLPNN Poly-SVM [CHAID MLPNN Poly-SV.. 

Figure 5. Stream of classification models for credit scoring system. 



if (AD in [ f ]) 

if (A13 in [ [JD]) 

if (A1 □ in [ 3D]) □□ Risky (cases 220; accuracy 0.945) 
if (AID in [ □□]) I Risky (cases 65; accuracy 0.985) 



if(A13in[ P}>q) 
if (A13 in[ BQ) 
if (AD in [ □□]) 

if (AID in [ 3D]) 
if (A15 



□□Creditworthy (cases 8; accuracy 0.625) 
□□Risky (cases 36; accuracy 0.861) 



□□i§ 

if (A14 .001 12 q 

if (A4 in [ OiD]) 
if (A4 in [ tyU]) 
if (A14 □ 12Dand A14 CD1DD) 
if (A14 □ 1 DDand A14 OD3E5) 
if (A13 in [ [JD] ) 
if (A13 in [ BID] ) 
if (A14 U3U5) 

if (A15 and A 1 5 □□ffid) 

if (A4 in [ OiD]) 

if (A14 □□ 1 5 CJ, 
if (A14 □ 15D) 
if (A4 in [ qrD]) 
if (Ai5 ftcrqi 



fjgj Creditworthy (cases 27; accuracy 0.926) 
□□Risky (cases 9; accuracy 0.667) 

□□Risky (cases 12; accuracy 0.833) 

□□Creditworthy (cases 21; accuracy 0.619) 
□□Risky (cases 10; accuracy 0.7) 

□□Risky (cases 12; accuracy 0.833) 



□□Risky (cases 7; accuracy 0.714) 
□□Creditworthy (cases 8; accuracy 0.EJ5) 
□□Risky (cases 8; accuracy 0.75) 
□□Creditworthy (cases 19; accuracy 0.947) 



if (AID in [ □□] ) 



if (A15 gp4) 

if (All ]jg4) 
if (All 04 ) 

if (A15 D4 and A15 #3237) 
if (A4 in [ OiD]) 
if (A4 in [ yD]) 
if (A15 D237 and A15 
if (A15 d2,ddq 



f ] Creditworthy (cases 32; accuracy 0.781) 
^Creditworthy (cases 28; accuracy 1.00) 



pj Creditworthy (cases 32; accuracy 0.781) 
CE Risky (cases 8; accuracy 0.5) 

□□2,1 1 I I) □□ Creditworthy (cases 84; accuracy 0.964) 
1 Creditworthy (cases 44; accuracy 1.00) 



Figure® The generated rule set of CHAID-DT model; this rule set sufficient for correct prediction of 89.88% of 

the training samples and DD7DD of the test ones. 



MLPNN classifier node is trained using the 
pruning method. It begins with a large network 
and removes the weakest neurons in the hidden 
and input layers as training proceeds. To 
prevent overtraining, portion of the training 
subset has been used as a validation subset. The 
network is trained on the rest of training subset, 
and accuracy is estimated based on the 
validation subset. The stopping criterion is set 



based on timeD the network is trained for one 
minute. However, the training process may be 
interrupted at any point to save the network 
model with the best accuracy achieved so far. 
Using the Australian credit scoring dataset, the 
resulting structure consists of four layers □ one 
input, one hidden layer and the output layer with 
3q 2 and 1 neurons respectively. The prediction 
accuracies of training and test samples are 
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□ .733| j and ipBDlp respectively. These results 
assured the conclusion of [2 J where they 
empirically stated that given the right size and 
the structure, the MLPNN is capable of learning 
arbitrarily complex nonlinear functions to 
arbitrary accuracy levels. 

Poly-SVM classifier node is used to train SVM with 
polynomial kernel. This model has four parameters 
that need to be optimized DC, J , r and d. The value of 
the regularization parameter C should be set between 
1 and 10 inclusive; increasing the value improves the 
prediction accuracy for the training data, but this can 



also lead to overfitting. Using trial and error we found 
that the best values for these parameters are IfJ 1, dl 
and 1 1 for C, y , r and d respectively. These values 
result in D5.D5D and ED41D prediction accuracies for 
training and test subsets respectively. It takes only 27 
seconds to build this SVM model. 

Filter, Analysis and Evaluation nodes are used to 
select and rename the classifier outputs in order to 
compute the performance statistical measures and to 
graph the evaluation charts. Table 2 summarizes the 
prediction results for the training and test subsets. 



Table 2. Prediction results for training and test subsets of Australian credit dataset. 



Model 






Dataset Partition 




Training 


Test 


Poly-SVM 


Correct 
■p rong 


□5.I3£P 

14.05a 


Q041 

13.5m 


MLPNN 


Correct 


mity. 


maij 


Jp rong 


mm 


mim 


CHAID-DT 


Correct 


mono 


m7m 


rong 


mi2ifc*i 


00220 




Figure 7. The cumulative gains charts of the three models for training and test datasets. 



The gains chart provides a visual summary of 
the usefulness of the information provided by 
the models for predicting categorical dependent 
attributes. The chart summarizes the utility that 
can be expected by using the respective 

predictive models, as compared to using baseline 
information only. The gains charts for Poly- 
SVM, MLPNN and CHAID-DT models trained 
in SPSS Clementine for training and test subsets 
are shown in Figure 7. The higher lines in the 
gains chart indicate better models, especially on 
the left side. These charts depict that the 
performances of the CHAID-DT and MLPNN 
models are similar for training and test samples 
and they are better than the performance of Poly- 
SVM. CHAID-DT model in the best in 

predicting the test samples with ED7EI 

prediction accuracy. MLPNN model is the 
second best model with HHHD accuracy. Poly- 
SVM model has achieved 71 .41 accuracy. 

Naturally, these percentages will increase if the 



number of training samples increases. □ hen the 
dataset is partitioned by the ratio LD2CD for 
training and test subsets, the CHAID-DT model 
achieved C2.57D prediction accuracy on test 
samples. 

Sensitivity analysis provides information about 
the relative importance of the input attributes in 
predicting the output attribute(s). The basic idea 
is that the inputs to the classifier are perturbed 
slightly, and the corresponding change in the 
output is reported as a percentage change in the 
output [23]. The first input is varied between its 
mean plus (or minus) a user-defined number of 
standard deviations, while all other inputs are 
fixed at their respective means. The classifier 
output is computed and recorded as the percent 
change above and below the mean of that output 
channel. This process is repeated for each and 
every input attribute. Figure □ shows sensitivity 
analysis of the three models in a graphical 
format. This graph states that four input 
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attributes have the most importance; A9, All!} 
A13 and A 14. However, AD is the most 
important one with a remarkable difference from 



the others. Also CHAID-DT model uses only six 
attributes out of the available fifteen to conduct a 
decision about the credibility of a consumer. 




5. Conclusions 

The objective of credit scoring is to help credit 
providers quantify and manage the financial risk 
involved in providing credit so that they can make 
better lending decisions quickly and more objectively. 
The decision is dependent on sociodemographic and 
financial characteristics of the consumer (such as age, 
education level, occupation and income), the 
repayment performance on previous loans and the 
type of loan. The automation of this process will 
result in an increase of the speed and consistency of 
the loan application process. □ ith the help of 
automated credit scoring, financial institutions are 
able to quantify the risks associated with granting 
credit to a particular applicant in a shorter time. This 
study evaluated the performance of three data mining 
methods □ support vector machine with polynomial 
kernel, multilayer perceptron neural network with 
pruning parameters and decision tree with Chi- 
squared automatic interaction in the prediction of 
credit scoring. The dataset contains information from 
real world credit bank. The dataset was partitioned by 
the ratio 7D3 Cgf for training and test the models. The 
CHAID-DT model achieved the best prediction 
accuracy with X7X success on the test samples. 
This accuracy increased to H2.57D, when the dataset 
is partitioned by the ratio □D2GP. CHAID-DT 
model is very fast and as all decision trees ones, can 
be expressed as a set of decision rules. It has the 
ability to pick and choose among the predictive 
attributes. Here, it used only six attributes out of 
fifteen to predict the output. These results recommend 
that CHAID decision tree model can be used 
effectively in real world systems to predict the 
credibility of loan consumers. 
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Abstract 

This paper shows the application of fuzzy logic 
based duty ratio control technique to reduce torque 
ripple in an induction motor employing Direct 
Torque Control (DTC). This technique increases the 
number of voltage vectors beyond the available 
eight discrete voltage vectors without any increase 
in the number of semiconductor switches in the 
inverter. In addition, this work incorporates Model 
Reference Adaptive System based observer to 
estimate rotor speed. Look-up table based on-line 
tuning PI controller is proposed for outer speed 
control loop to achieve swift response, less 
overshooting and precision speed control to have 
wide torque-speed characteristics. A new algorithm 
for optimized value of stator flux based on the 
maximum reference value of electromagnetic torque 
is proposed to operate in conjunction with duty 
ratio control. The performance of the proposed 
drive system is evaluated through digital simulation 
using MATLAB-SIMULINK package. The simulation 
results clearly depict the superiority of devised 
method over the existing methods of DTC. 

Keywords: Direct Torque Control , Fuzzy Logic Duty 
Ratio Control , Induction Motor , Model Reference 
Adaptive System. 
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Instantaneous values of rotor flux 
linkages DQ-axes components in the 
stationary reference frame, (wb) 
Instantaneous values of stator flux 
linkages DQ-axes components in the 
stationary reference frame, (wb) 

Stator and rotor resistance, Q 
Stato, rotor and mutual inductance, H 
Electromagnetic and load torque, N-m 
Reference electromagnetic torque, N-m 

Maximum value of reference 
electromagnetic torque, N-m 
Number of poles 
Inductance leakage factor 
Duty ratio 

Acceleration, radian/sec 2 
Degree of membership 
Slip speed, radian/second 
Stator flux, Wb 

Stator flux and current space vector 
Absolute value of Stator flux 
space vector, Wb 
Rotor transient time constant 

Torque and flux error 



1. Nomenclature 

v„,v„ Instantaneous values of the rotor 

voltage DQ- axes components in the 
stationary reference frame. (V) 
v„,v„ Instantaneous values of the stator 

sD ’ sQ 

voltage DQ- axes components in 
the stationary reference frame. (V) 



2. Introduction 

The Induction Motor (IM) sensorless speed control 
replaces the speed sensors by microcomputer or 
DSP based estimation, which uses machine terminal 
voltage and current, and thus improves the 
robustness and reliability of the drive with reduction 
in cost and complexity. The Model Reference 
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Adaptive System (MRAS) based speed estimation 
method is proved as one of the best due to its high 
performance ability and straightforward stability 
approach [25]. The IM drives controlled with the 
vector control method has found wide acceptance 
in the industry. However, this control technique 
requires complex coordinate transformation, inner 
current control loop and accurate system 
parameters [3]. On the other hand, the Direct 
Torque Control (DTC) method [7, 12] provides 
robust and fast torque response without such 
coordinate transformation, PWM pulse generation 
and current regulators [20]. Moreover, DTC 
minimizes the use of motor parameters [7, 10]. 
Despite, this technique suffers from a major 
disadvantage of steady state ripple in torque and 
flux, because none of the inverter-switching vector 
is able to generate the exact stator voltage at 
proper instants as well as in space. These torque and 
flux ripples affect the accuracy of speed estimation; 
result in high acoustic noise and harmonic losses 
[26]. There are many methods to reduce this torque 
and flux ripple: (a) the alternative inverter topologies 
[11], multilevel inverters [6, 28] and matrix 

converters [4] which increase the number of 
switches, and thus cost and complexity; (b) the 
higher switching frequencies reduce the harmonic 
content of stator current and thus torque and flux 
ripple. However, such higher switching frequencies 
lead to increased switching losses and stress on 
semiconductor switches of the inverter [22, 23]; (c) 
yet, another method of reducing torque and flux 
ripples is space vector modulation [9, 16, 18] that 
has disadvantage of variable switching frequency 
[14]. Moreover this method requires complex 
control schemes than classical DTC and is machine 
parameter dependent; (d) the Discrete Space Vector 
Modulation [5, 10] overcomes the disadvantages of 
SVM technique with an accurate switching table and 
five level hysteresis bands, but this method cannot 
guarantee its functionality at low speed range, 
especially with heavy load [20]. 

Some more solutions available in literature include 
modified switching table based DTC [15], variable 
amplitude control of flux and torque hysteresis 
bands [13], open loop control of hysteresis band 
amplitude [27], and fuzzy logic based variable 
amplitude control of flux and torque hysteresis 
bands [8]. But these are complex in nature for easy 
implementation. The ripple in the torque and flux 
can be easily reduced by applying the selected 
inverter vector only for the part and not for the 



entire switching period unlike that in the classical 
DTC. This technique, also known as duty ratio 
control, increases the number of voltage vectors 
beyond the available eight discrete ones, without 
any increase in the number of semiconductor 
switches in the inverter [21, 23]. The work in [2] 
refers to fuzzy logic duty ratio control but fails to 
give any details. Similar work presented in [19, 24] 
claims reduction in ripple to one-third value. 
However, simulation results were shown for low 
speed operation only using a fuzzy controller with 
three membership functions for inputs and output. 
They have shown an algorithm for optimized value 
of reference flux based on reference load torque. 

The work presented in this paper considers the 
design of new fuzzy logic controller with five 
membership functions for the inputs & output 
along with the optimized flux algorithm based on 
maximum reference value of the electromagnetic 
torque to adjust the "duty ratio" of inverter 
switching vectors. The present work shows much 
better performance (reduced ripple) though it has 
some similarities with those in [2, 19, 23]. A 
significant reduction in the ripple than that in [19, 
24] has been achieved in the present work, 
incorporating MRAS and adaptive PI controller. The 
sensorless mode of operation has been achieved 
through MRAS based speed estimation. The 
adaptive PI controller that is based on look-up table 
is used for outer speed loop for precision speed 
tracking. A series of simulation tests are conducted 
using MATLAB-SIMULINK package to validate the 
performance of the devised algorithm. 

The complete paper is organized as follows: Section 

3 explains the modelling of induction motor, Section 

4 explains strategy of torque ripple minimization. 
Section 5 discusses design of fuzzy logic duty ratio 
controller. The simulation results, comparison and 
discussion are presented in Section 6. Section7 
concludes the work. 

3. Modelling of Induction Motor 

In DTC the induction motor is modeled in stationary 
reference frame, thus it does not require any 
transcendental equations as they require in FOC. 
Stator and rotor fluxes can be expressed as follows. 

Ad = J ( V sD ~ R s To \dt 

A'sQ=\( V sQ- R s- i s Q \ dt 

Ad = J { V rD - KId ~ P M mA<Q \ dt = J (“A Ad “ P ^mA Q )- dt 

Ae = j [ V rQ ~ R rl Q + P-^mAo \ dt = j {~ R r K Q + P ^mAo)- dt 

[For squirrel cage induction motor v rD = v rQ =0] 
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Stator and rotor flux linkage equation can be written 
as follows: 



A- <2 AAg + L m i r Q 

^ sD ~ LJsD ~l~ AwA) 

Kq — AAg + AAq 

Ad ~ AAd An A) 



Solving the upper set of equations the stator and 
rotor current equations are derived as follows: 
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whereL = L .L - L“ 



(3) 



reduced by applying the selected inverter vector 
only for the part and not for the entire switching 
period unlike that in the classical DTC IM drive. The 
time for which an active voltage vector has to be 
applied is chosen just to increase the 
electromagnetic torque to its reference value and 
the null voltage vectors are applied for the rest of 
the switching period. During the application of the 
null vectors no power is absorbed by the motor, and 
thus the electromagnetic torque is almost constant 
or decreases slightly. But this decrease should be 
small to have minimum torque ripple. Since this 
decrease in torque also depends on the modulus of 
the reference stator flux, an optimized value for the 
flux has to be used, which is large enough to 
generate the reference torque. This implies that the 
maximum electromagnetic torque reference has to 
be found and the optimized stator flux reference 
corresponds to this, as follows: The expression of 
stator current space vector in stator flux oriented 
reference frame can be obtained as, 

r : (L m /L s )\jT s \(R r + jca sl L r ) 

{R r + jco sl oL r ) 



Electromagnetic torque expression of machine is as 
follows: 




The speed can be calculated from the following: 

T e ~T l = J-~~ + B-co m (5) 

at 



4. Torque Ripple Minimization Strategy 

Despite of the available complex solutions [4-6, 9- 
11, 14, 16, 18, 26, 28] the duty ratio scheme presents 
best remedy to minimal torque and flux ripple, by 
overshadowing the above-mentioned existing 
methods. In the classical DTC, a voltage vector is 
applied for the entire switching period, and this 
causes the stator current and electromagnetic 
torque to increase over the whole switching period. 
Thus for small errors, the electromagnetic torque 
exceeds its reference value early during the 
switching period, and continues to increase, causing 
a high torque ripple. This is then followed by 
switching cycles, where the null switching vectors 
are applied to set the electromagnetic torque to its 
reference value. 

The ripple in the torque and flux can be easily 



The imaginary part of the above equation gives the 
torque producing stator current ) , and maximum 

value of this occurs at tD sl max = 1/L r ' (T' denotes the 

rotor transient time constant). By substituting the 
corresponding i sy mSiX value in the electromagnetic 

torque expression, T e = 3/2(P/2)y/ s xi s of induction 

motor, the maximum torque can be evaluated as 
follows: 

T e,M a = 3 / 2 { P / 2 )\V's\i S y,rmx (?) 

The maximum electromagnetic reference torque can 
be found from the above equation by 
substituting T e = T e ref , as shown in Equation (3). 
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The optimized reference flux linkage for the given 
maximum electromagnetic torque reference can be 
calculated by setting \y s , re f = \y/s as ' n Equation (4). 
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( 9 ) 



When the above optimized stator flux linkage value 
is used in DTC induction motor drive along with 
duty ratio control, the torque ripples reduce 
significantly. But, the duty ratio of each switching 
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period is a non-linear function of the 
electromagnetic torque error ( E Te ), stator flux- 
linkage error ( E l//S ), and the position of the stator 
flux-linkage space vector (p s ). It is difficult to model 
such non-linear function. However, the 
characteristics (such as, model free nature and non- 
dependence on mathematical equations [17]) of 
fuzzy logic controller makes the duty ratio 
determination possible and easier during every 
switching period. 

5. Proposed Methodology 

Fig 1 shows the schematic of proposed DTC IM 
drive with fuzzy logic based duty ratio control. In 
this paper, two Mamdani type fuzzy logic controllers 
(FLC-1 and FLC-2) that contain fuzzifier, inference 
engine, rule base, and defuzzifier are chosen. The 
FLC-1 and 2 are used for positive and negative flux 
error respectively. In classical DTC, outer loop speed 
regulators are conventional PI controllers, which 
require precise mathematical model of the system 
and appropriate value of PI constants to achieve 
high performance drive. Therefore, unexpected 
change in load conditions or environmental factors 
would produce overshoot, oscillation of the motor 
speed, oscillation of the torque, long settling time 
and, thus causes deterioration of drive performance. 
To overcome this problem, a look-up table is 
designed from the experiences of speed response of 
classical DTC. According to the speed error and 
change of speed error, the proportional and integral 
gains are adjusted on-line. The speed can be 
calculated by the MRAS [1, 25], which uses the 
motor terminal voltage and current as shown in Fig 
1. 




Figure: 1. Schematic diagram of proposed DTC with fuzzy logic 
duty ratio controller 



A. Selection of Input/output Variables 

The design starts with assigning mapped 
input/output variables of FLC. In this work, the first 
input variable is the torque error ( E Te - T e re f - T e ) 

and the second input variable is the stator flux 
vector position ( p s ) at a sampling tinner^. The 
output variable is duty ratio {S). 

B. Fuzzification 

The success of robust fuzzy controller design mainly 
depends upon this stage. In this stage the crisp 
variables of the inputs E Te (t s ) and p s (t s ) are 

converted into fuzzy variables E Te and that can 
be identified by the levels of membership in the 
fuzzy set. Each fuzzy variable is a member of the 
subsets with a degree of membership p varying 
between 0 (non-member) and l(full member). 

To make the torque and duty ratio variations 
smaller, the universe of discourse of torque error 
and duty ratio are divided into five overlapping 
fuzzy sets. However, to reduce the complexity of 
design, the stator flux position is defined with three 
overlapping fuzzy sets only. The universe of 
discourse of all the variables, covering the whole 
region is expressed in per unit values. The fuzzy 
subsets are defined with triangular membership 
functions as shown in Fig 2. The linguistic labels are 
defined as VS=Very Small, S=Small, M = Medium and 
L=Large, VL=Very Large. 




(a) (b) 



Figure 2: Membership functions for (a) torque error & duty ratio 
(b) stator flux position. 

C. Rule Base and Fuzzy Inference Engine 

The Rule base involves defining the rules in IF-THEN 
form, which govern the relation between input and 
output variables in terms of membership functions. 
In this stage the input variables are processed by 
the inference engine that executes 15 rules (5x3) as 
shown in Table 1. These rules are generated from 
the knowledge of control systems and the 
simulation results of classical DTC induction motor 
using different switching states. 
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The inference engine also includes the application 
of fuzzy operator (AND, OR), product operation of 
fuzzy implication and maximum aggregation. The 
relation between different conditions in the same 
rule is done by means of 'AND' operator. On the 
other hand, the relationship between different rules 
is done by means of 'OR' operator. 



Table 1: Fuzzy control rules for duty ratio determination 
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D. Defuzzification 

This stage introduces different inference methods 
that can be used to produce the fuzzy set values for 
the output fuzzy variables. In this paper, the center 
of gravity (COA) or centroid method is used to 
calculate the final fuzzy value. The COA expression 
with a discretized universe of discourse can be 
written as, 

Unit (A ) 

S = ~n ( 5 ) 

'y ' i ^out ( 4 " ) 

;=i 

Where, n denotes the total number of rules. The 
comparator circuit as shown in Fig 1 compares the 
duty ratio determined during each switching period 
with a triangular signal, whose period is equal to 
that of switching period and thus determines the 
duration for which active vector should be applied. 
The modified symmetrical switching vectors fed to 
the inverter would improve the performance of 
drive. 

6. Results, Comparison and Discussion 

A series of simulation tests are conducted on a 4KW, 
4pole sensorless inverter-fed IM to evaluate the 
performance of proposed DTC method. 

The simulation results for speed response as well as 
speed estimation by MRAS are analyzed during 
transient and steady state conditions. However, the 
torque response has been analyzed in steady state. 



These tests are carried out for high (220rad/sec) and 
low (50rad/sec) speed operation with sudden drop 
in load torque from 20N-m to 5N-m at 1.5sec. A 
constant un-optimized stator flux value of 0.5 Wb is 
used for classical DTC. Whereas, optimized stator 
flux values of 0.3734Wb (for 20N-m load) and 
0.1867 Wb (for 5N-m load) are used for proposed 
one. 





Time, sec 



Time, sec 



(c) 



(d) 



Figure 3: Simulation results for sensorless speed estimation and 
response due to (a) classical DTC with reference speed of 
220rad/sec (b) proposed DTC with reference speed of 220rad/sec 
(c) classical DTC with reference speed of 50rad/sec, and (d) 
proposed DTC with reference speed of 50rad/sec. 




Figure 4: Simulation results for electric torque response due to (a) 
classical DTC with reference speed of 220rad/sec (b) proposed 
DTC with reference speed of 220rad/sec (c) classical DTC with 
reference speed of 50rad/sec, and (d) proposed DTC with 
reference speed of 50rad/sec 

Figure 3(a) and 3(b) shows reference speed of 
220rad/sec, measured speed, estimated speed, error 
between measured and estimated speed, and load 
torque on a single diagram for classical DTC and 
proposed DTC respectively. Figure 3(c) and 3(d) 
show similar result for a reference speed of 
50rad/sec. The smoother response for 3(b) and 3(d) 
against 3(a) and 3(c) respectively, reflects ripple 
reduction in torque, flux and stator current. As 
shown in figures, the speed estimator is able to 



25 



ACSE Journal, ISSN: 1687-481 1 , Volume 10, Issue 1, ICGST LLC, Delaware, USA, November 2010 



track the rotor speed under all-operating 
conditions. The smoother and quicker response with 
the on-line tuned PI controller for all speed-torque 
operating range validates its functionality. 

Figure 4(a) and 4(b) shows electric torque response 
for classical DTC and proposed DTC respectively for 
high speed operation. Figure 4(c) and 4(d) show 
similar results for low speed operation. In addition 
to the inherent disadvantage of classical DTC, the 
constant reference flux causes higher ripple at lower 
torque level because, the chosen 0.5 Wb is larger 
than the optimized flux value (0.1867Wbs). But, the 
proposed method causes lower ripple at lower 
torque level than at the higher torque level. 
Moreover, this method is able to eliminate the 
torque undershoot and overshoots. The average 
amount of ripple reduction in torque compared to 
classical DTC is observed to be l/19 th (~llN-m for 
classical with peaks up to ~16N-m, and ~1.2N-m for 
proposed) and l/15 th (~15N-m for classical with 
peaks up to ~19N-m, and ~lN-m for proposed) 
respectively for higher and lower torque levels of 
high-speed operation. Similarly, for low speed 
operation the reduction is observed to be more 
than l/10 th (~9N-m for classical with peaks up to 
~14N-m, and ~0.8N-m for proposed) and l/17 th 
(~lll\l-m for classical with peaks up to ~15N-m, and 
~0.6N-m for proposed) respectively for higher and 
lower torque levels. 

7. Conclusion 

In this paper, three current research topics viz. 
sensorless control, on-line PI tuning and torque 
ripple minimization strategy are brought to a 
common platform to alleviate the disadvantages of 
classical DTC employed to control an inverter-fed 
induction motor drive. The fuzzy logic duty ratio 
controller with optimized reference flux command is 
applied and verified. Low speed operation of the 
proposed method dictates the proper functionality 
of speed estimation as well as precision torque 
control. From the simulation tests the following 
justifiable conclusions are drawn against the existing 
solutions. 

1. The devised method is easier to understand, 
design and implement. The input/output scaling 
factors for fuzzy logic controller are absent. 

2. Increased efficiency as well as lower acoustic 
noise for the induction motor drive is observed 
because of the significant reduction (~l/15 th ) in 
torque, flux and current ripples/harmonics. In 
addition to the torque ripple minimization, under 
(over) shoots are also eliminated. 



3. The proposed method operates at a lower 
switching frequency, and thus reduces switching 
losses as well as stress on semiconductor 
switches of the inverter. 

4. The propsed algorithm is able to track the actual 
speed under all operating conditions. 

5. Durig transient conditions the error between 
measured and estimated speed for fuzzy MRAS 
is observed to be low when compared to MRAS. 



Appendix 

3-Phase Induction Motor Parameters 

Rotor type: Squirrel cage, Reference frame: 

Stationary 

4KW, 1440 rpm, 50Hz, 4 Poles, R s = 1.57Q, R r = 1.21 

Q, L s = 0.17 H, L r = 0.17 H, L m = 0.165 H, J = 0.06 

Kg-m 2 .. 
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Abstract 

During our studies of the techniques of conception 
(design) and modelling of linear and nonlinear 
microwave circuits, we noticed the big need to a 
method able to improve the analysis and the 
understanding of circuits, which often function in 
high frequency. Therefore, we propose, in this paper, 
a new methodology to study a microwave filter with 
localized elements (Chebychev low-pass filter) by a 
joint application of the bond graph approach and 
the scattering formalism. This study consists of, on 
the one hand, determining and simulating the 
scattering parameters of this filter and on the other 
hand, to modeling the incident and reflected wave 
propagation since the source towards the load 
through this filter. In fact, the modelling of these 
various wave propagation amounts modelling the 
scattering matrix constituted by the scattering 
parameters of the studied filter on a particular type 
of bond graph model often named: "Scattering 
Bond Graph". 

Keywords: Bond graph modelling, scattering 

formalism, Chebychev low-pass filter, power wave 
propagation, and simulation 

1. Introduction 

The scattering formalism [1] was used in vast 
physics fields such as the characterization of the 
electric circuits. Several work, since the invention of 
the bond graph approach [2], showed that the 
scattering formalism constitutes an alternative 



approach for the physical systems modelling [3]. 
They pointed out some properties, and in particular 
the orthogonality of the scattering matrix [4] or 
wave matrix [5] which respects intrinsically the 
causal relations and includes explicitly the 

conservation laws [6]. In addition, they showed that 
the scattering representation exists for systems 
having neither impedance nor admittance such as 
the junctions of Kirchhoff, the gyrateurs and the 
transformers [7]. 

Moreover, the bond graph language [2] is based on 
a graphic representation of the physical systems. 
These representations are based on the 
identification and the idealization of the intrinsic 
characteristics of the physical environments and on 
the structuring of a complex physical system in the 
networks form [8]. Besides, in physics, the analogies 
theory allows bond graph technique [2] to 
generalize the representation networks with all the 
traditional physics fields of the systems with 
localized and/or distributed parameters [9]. 

The proposal of this paper is to apply the scattering 
formalism [1] and the bond graph technique [2] 
jointly to build a model able to highlight the 
incident waves and reflected waves and their 
propagation on a particular type of bond graph 
model often named "scattering bond graph" of a 
Chebychev low-pass filter which functions in high 
frequency. 

At first, we propose to determine the scattering 
parameters (transmission and reflexion coefficients) 
of the studied system from its transformed and 
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reduced bond graph model [10] by the integro- 
differentials operators [11] which based on the 
causal ways and algebraic loops present in the 
causal bond graph model Then, we will make a 
simple comparison by simulation of the found 
scattering parameters and the classic techniques of 
conception and simulation of the microwave circuits 

[12] at the aim to validate the founding results. 
Finally, we will apply to these scattering parameters 
a procedure developed and described thereafter to 
obtain the famous "Scattering bond graph" model 

[13] . 

2. Study of LC filters with localized 
elements 

The LC filters are passive filters including only 
inductances and condensers, in the quadripole form 
inserted between two terminal resistances: the 
generator resistance Ri, and the load resistance R 2 . 
There are many topologies to produce these filters 
of which the most answered is that the LC filters in 
scale. The LC filters in scale are not any more to 
show, which is based on the knowledge of the 
transfer functions and characteristic function 
corresponding to a specified template [14]. 

Starting from a transfer function and a characteristic 
function corresponding to a specified template, it is 
enough to calculate the various elements L and C 
for our structure, calculate the entry impedances 
seen from the two ends of the filter, then, extract 
the elements (with standardized (normalized) value) 
in an iterative way. The elements values are 
calculated by the method known as Darlington 
method [15]. 

2.1 Synthesis of LC filter with localized elements 

We propose to synthesize a low-pass filter with 100 
MHz of cut-off frequency (f c ), sensibility: k = 0.5, and 
attenuation: A max = O.ldB, A min = 20dB. 

By choosing a Tchebytchev approximation function, 
we obtain a filter with order 4 whose transfer 
function is as follows: 

H (P)= 4 y 1 2 (1) 

1.206/7 +2.177 p +3.170/? +2.445 p+l 

It is about an asymmetrical structure with even 
order. The extraction according to the Darlington 
method [15] leads to the standardized (normalized) 
values below, and corresponds to a n (pi) structure 
as figure 1 refers to. 
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Figure 1: The corresponding Chebychev low-pass filter. 



2.2 Simulation results of the scattering 
parameters 

The traditional tool for simulation of the scattering 
parameters under HP-ADS software [12], as figure 2 
refer to, makes it possible to visualize the filter 
response in transmission (Si 2 and S 2 i) and in 
reflexion (Sn and S 22 ) on a band frequency 
surrounding the cut-off f c as figure 3 refers to. 
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Figure 2: Tchebytchev low-pass filter, under HP-ADS software 
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Figure 3: Simulation results of the Tchebytchev low-pass filter. 



3. Scattering parameters of the filter 
exploited from its bond graph model 



We propose in this part to determine the scattering 
parameters [4] of the Chebychev [16] low-pass filter 
from its transformed and reduced bond graph 
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model [10] and by paying attention to the causality 
assignment and to compare the obtained results 
with the results found previously by the traditional 
techniques used for the conception of the linear and 
nonlinear microwaves circuits under HP-ADS 
software [12]. 

3.1 Relations between wave-scattering matrix 
and bond graph model 

Generally, any physical system exists in the form of a 
quadripole inserted between two particular ports Pi 
and P 2 which respectively represent the entry 
(source) and the exit (load) of the total system [17]. 
This system can be represented by a generalized 
bond graph model transformed and reduced as the 
figure 4 indicates it. 



Port 




Port 



Process i 
/ 



Port 2 



Port 



Figure 5: Different possibilities of causalities assignment. 

For each type of reduced and causal bond graph 
model given by the above figures, we will 
respectively have the following matrices. 



<P 2 



(«U h 12 



Figure 4: General transformed and reaucea oond graph model. 
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And to establish the entry-exit analytical relations, 
the bond graph model of the studied system must 
be transformed, reduced and especially be causal 
since these relations rest on the concepts of causal 
way and algebraic loops which exist in the reduced 
bond graph model [10, 17], 

The causality assignment to the reduced bond 
graph model of figure 4 enables us to notice that 
they are four different cases of causality assignment 
in input-output of the process [17], 



(3) 

(4) 

(5) 

( 6 ) 



From the relations given above, we can define the 
matrix often noted H which gathers the various 
integro-differentials operators Hy. 

H U #12 



#= 



h 2 \ H 22 



(7) 



Hy represent the integro-differentials operators 
associated to the causal ways connecting the port Pj 
to the port Pi and obtained by the general form 
given below. 

£ GkAk 
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A=l -I ^ +y LjLj -y LjLj L k +■ ■ -+(-l) m Z- • • 

Where: 

• A = the determinant of the causal bond graph. 

• Hjj= complete gain between Pj and Pj. 

• Pi = input port. 

• Pj= output port. 

• N = total number of forward paths between Pj and 



' —A 



Port 2 



Pj- 



• Gk = gain of the k th forward path between Pj and Pj. 

• Lj = loop gain of each causal algebraic loop in the 
bond graph model. 



cp 2 



Case (b) 
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• LjLj = product of the loop gains of any two non- 
touching loops (no common causal bond). 

• LjLjL k = product of the loop gains of any three 
pairwise nontouching loops. 

• A k = the cofactor value of A for the k th forward 
path, with the loops touching the k th forward path 
removed; i.e. Remove those parts of the causal bond 
graph which form the loop, while retaining the parts 
needed for the forward path. 

We noted that: 



a ,=f^L , b . = S J^L 

l 2 i 2 



£. = 
l 
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( 10 ) 
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These are the reduced voltage and current with 
respect to R 0 (scaling resistance). 
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According to the preceding equations, we can find, 
for each case of causality assignment which 
indicated by figure 5, one wave matrix. 




(15) 



(16) 



(17) 



(18) 



AH=Hnff 2 2-tf 12^21 (19) 

The scattering parameters are given by the 
following scattering matrix: 



=[ 5 ] 



( 20 ) 



The relations between these matrixes are given by 
these following equations: 



w i2= s n 



-l 



W 2\~ S \2 S ll* S 22* S 21 
W 22 =S ll* S 2l 



( 21 ) 



And the corresponding scattering matrix is given 
below: 



5 (r ) = 



W 22 * W 1 2 W l 1 _W 2 1 * W 1 2 * W 22 



Wr 



-1 



22 



~ W 21* W 22 



( 22 ) 



3.2 Application on the Chebychev low-pass filter 

Here, we will try to apply and check the procedure 
described previously on the Chebychev low-pass 
filter [16] given by figure 1. 

The bond graph model of this filter is given by 
figure 6. 




20*im4. 1 Viewer (c) CLP 2009 



Figure 6: The corresponding bond graph model. 

To extract the scattering parameters [4] from the 
bond graph representation and by using the new 
method which is described previously; we must 
transform the bond graph model given by figure 6 
into a causal bond graph model often named 
reduced bond graph model [10] only containing the 
decomposition junction (1-ionction or 0-ionction) 
and the reduced variables with respect to a scaling 
resistance R 0 (internal source resistance). 



Decomposition 




Figure 7: The reduced bond graph model. 
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Where: 
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By decomposition the reduced bond graph given by 


figure 7, 


we will have the following bond graph 


representation: 
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Figure 8: The first sub-model 


Figure 9: The second sub-model 


Notice that: 






• z\ : the 


reduced equivalent impedance of the i 


element put in series. 






• yj: the 


1 reduced equivalent admittance of the i 


element put in parallel. 






So we have: 
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p: The Laplace operator. 

The tow sub-models are in conformity with case (a) 
described previously. So we have the integro- 
differentials operators by taking into account to the 
previously equations. 

1 ^=— : Loop gain of the algebraic loop given by 

z i y i 

the first sub-model. 

l 2 =^L : Loop gain of the algebraic loop given by 

* 2^2 

the second sub-model. 
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the second sub-model. 



From these operators, we can deduce directly the 
wave matrix of the first and second sub-model by 
taking into account to equations of case (a) : 



2 

2 



z \ y \~ z \~ y \ +2 ~ z l y l +z l +y l 
~ z l y l~ z l +y l z \ y \ +z \ +y \ 

z 2 y 2~ z 2~ y 2 +1 ~ z 2 y 2 +z 2 +y 2 
~z 2 y 2 ~z 2 +y 2 z 2 y 2 +z 2 +y 2 



(29) 

(30) 



The wave matrix of the complete system can be 
given by the product of the first and the second 
wave matrix such us: 



=wO) *w^ = 



(31) 



From this matrix we can deduce these following 
scattering parameters: 
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A, =1+ — ( — : Determinant of causal bond graph of the 
first sub-model. 



A 2 =l + ^ : Determinant of causal bond graph of the 



z 2 y 2 

second sub-model. 



3.3 Simulation results of the scattering 
parameters 

A simple programming and simulation of the 
following scattering parameters equations, give the 
figure 10, figure 13, figure 11 and figure 12 below 
which represent respectively the reflexion and 
transmission coefficients of the studied filter. 
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5 n 



T c\ T c2 T L\ T LlP A+T Ll T c2 (T c\- T L2 ) P i + 
\_ T c\ ( ' T L2 +T Ll )+T c2 ( ' T L2~ T L0\ p2+( ' T c\ +T c2~ T Ll~ T L2 )p 

T cfc2 T Ll T L2 p4+T Ll T c2^cl +T L2^ p3 + 

^c\ +T c2 ) ^L\ + z L2 )p2+( ~ T c\ + T c2 +T Ll +T L2 )p+2 



(36) 



T cl T c2 T L\ T L2 p4 +T Ll T c2 ( ' T cl r L2 )p3 + 

[ r c l< r £2- r £l>- r C 2< r £2 +r Ll>]p 2+ < r cl +r C 2- r £r r £2 ) P 

T cl T c2 T Ll T L2 p4 + T Ll T c2 ( - T c\ +T L2' )p3 + 
( T cl +T c2'> ( - T Ll +T L2'> p2+( - T cl +T c2 +T Ll +T L2') p+2 



(37) 



2 

T c\ T c2 T L\ T L2 p4+T L\ T c2 i - T c\ +T L2'>P i + 
( ' T cl +T c2 )( ' T Ll +T L2 )p2+( ' T c\ +T c2 +T Ll +T L2' ,p+1 

2 

T cl T c2 T Ll T L2 p4+T Ll T c2 ( - T cl +T L2' lp3 + 
( - T c\ +T c2 )( - T Ll +T L2 )p2+< ' T cl +T c2 +T L\ +T L2' ,p+2 



(38) 



(39) 




Figure 10: Reflexion coefficient Sn seen at entry 




Figure 11: Transmission coefficient S i2 seen from exit to entry 



20 40 60 



100 120 140 160 

freq,MHz 



200 220 240 260 280 300 



Figure 12: Transmission coefficient S 2 i seen from entry to exit 




freq,MHz 



Figure 13: Reflexion coefficient S 2 2 seen at exit 

The purpose of these simulations is to validate this 
new determination method of scattering parameters 
from a causal bond graph model of a filter often 
functioning in high frequency and contrary to the 
work carried out by Pr A. KAMEL [3] where the 
causality concept, which represents a significant 
property in the formalism with network type and in 
particular the bond graph formalism, was ignored. 

All simulation given by figure 10, 11, 12 and 13 were 
carried out in the MAPLE software. They give us the 
curves of the reflexion coefficients (Sn and S22) and 
transmission coefficients (S i2 and S 2 i). These curves 
give us information about the cut-off frequency and 
the type of the studied filter which is a low-pass 
filter with cut-off frequency 100 MHz. 

The simulation given by the figure 3 which was 
carried out in the HP-ADS software and the 
simulation by the new analytical method given by 
figure 10, 1, 12 and 13 all show that the two 
representations are equivalent. Thus we can say that 
this new method is validated. It gives us always the 
same results given by the classic techniques of 
conception and simulation of the microwave circuits. 

4. Scattering Bond Graph realization 

The scattering matrix [4] of our studied process is 2- 
2. It is a particular form [18] despite the complexity 
of expression of the series impedance or parallel 
admittance [3]. It is orthogonal since the process is 
considered without loss, is not a transfer matrix and 
it admits the following general form. 



s= 



b h lpn + ‘ ' ‘ +b 0 bl n pU + ‘ ' ‘ +b 0 2 
b l l P n+ - +b l X b l 2 P n+ - +b n 2 . 



a p n +a A p n * +---+a n 
n y n-\ l 0 



We note that: 

d(p)=a n p n +a n _ x p n ~ l +■ ■ • +a 0 



(40) 



(41) 
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Now, we will use the method which is developed by 
Pr. A. Kamel. 

Indeed, if we consider the scattering matrix form 
found above for the process alone, we can say that 
it is not a true transfer matrix [8]. Moreover, it is not 
in the adequate form since its various Sjj parameters 
and sometimes Sy have the numerator's degree 
equal to that of denominator and that poses a 
major problem to determinate the scattering bond 
graph model of any physical system studied in a 
general way. 

The solution to this problem is to consider the 
scattering matrix of a process as a transfer matrix 
connecting the incident and the reflected waves in a 
symbolic system form. 

We start by carrying out an Euclidean division of 
each term of the numerator matrix (scattering 
parameters) by the common denominator d(p) 
which leads to the new shape of the scattering 
matrix [19] such as: 



s=s'+d (42) 

S: the new scattering matrix with degrees in the 
numerator at most one less than that of d(s). 

D: direct transmission matrix (The rest of the 
Euclidian division). 



elements whose modules are precisely the 
coefficients. 

• Add the direct part (transmission matrix D) by 
using information bonds. 

• To obtain the scattering bond graph model of 
the physical system, it is enough to add the 

reflexion coefficient P =^of the source and 

* ^o +1 

the reflexion coefficient Pn = of the load to 

C "L +1 

the scattering bond graph model of the 
process to the adequate sites. 

It is interesting to notice that the structure of the 
scattering bond graph of the process remains the 
same whatever the degree of the common 
denominator of scattering matrix. The only thing 
that changes is the corresponding number of I and 
C linked to the a-Routh expansion [3, 13]. 



2T c2 T Ll T L2 p3+2T Ll T c2 p2+2 ( T Ll +T L2^ p+2 
T cl T c2 T Ll T L2 p4+T Ll T c2^ T cl +T L2^ p3 + 
^cl +T c2^ T Ll +T L2^P 2+( ~ T cl +T c2 +T Ll +T L2^P +2 
2 

T cl T c2 T Ll T L2 p4+T Ll T c2 ( - T cl +T L2' >p3 + 
( - r c\ +T c2^ r Ll +T L2^P 2+< ~ r cl +T c2 +T Ll +T L2^P +2 

2 

T cl T c2 T Ll T L2 p4+T Ll T c2^cl +T L2'> p3 + 
(T cl +T c2^Ll +T L2^P 2+(T cl +T c2 +T Ll +T L2^P +2 



(44) 

(45) 

(46) 



d 4 ^2 



(43) 



Thereafter, we seek to apply the alpha-beta method 
of development in continuous fraction starting from 
the Routh method [20] and building the 
corresponding bond graph model, since it is about a 
multivariable system [18] while being based on the 
following systematic procedure: 

• Calculate the cx-Routh table from the common 
denominator d(p) and the (3-Routh table from 
the new numerator of the S -matrix. 

• Construct the direct chain by using the 
adequate number of elements I-C (which otj 
coefficients are their modules) in integral 
causality, equal to the degree of d(p). 

• Duplicate this chain and construct the two 
entries of the quadruple. 

• Construct the tow outputs by using information 
bonds and a sufficient number of TF and GY 



D= 



2T c\ T c2 T L\ p3+2T L\ T c2 pl+2 ^ T cl +T c2 )p+2 
T cl T c2 T Ll T L2 p4 +T Ll T c2 (r cl +r L2 )p?> + 
( T cl +T c2^ T Ll +T L2^ P 2 + ( T cl +T c2 +T Ll +T L2^ p+2 



1 0 

0 -1 



(47) 



(48) 



The alpha-beta coefficients are then: 
From d(p) we have: 

• ai, a 2l a 3 and a 4 . 

From S n we have: 

• A“- 02 > 03 and 04 



From S n = S 21 we have: 



A 12 = A 21 - 02 = 00 . A 12 = 00 and 



?21 



21 



i21 



00 = 00 



From S 22 we have: 



00 , A 22 - A 22 and A 2 
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Process 



20sim4.1 Viewer (c) OP 2009 

Figure 14: Scattering bond graph model of the low- 
pass filter connecting to its source and load 

5. Conclusion 

In this study, we showed how we can find the 
scattering parameters (scattering matrix S) of 
Tchebytchev low-pass filter functioning in high 
frequency by a jointly use of the scattering 
formalism and the bond graph technique and by 
taking into account the causality concept. 

Indeed, by adopting the reticulation assumption 
which allows the phenomena separation and the 
orthogonality property of the scattering matrix S 
associated to a physical system, we showed the 
existence of an ideal structure junction of Kirchhoff 
which includes a junction with common effort and 
flow known as a decomposition junction (0-junction 
or 1-junction: according to the studied filter type). 

The new determination method which applied to a 
filter functioning in high frequency, makes it 
possible, on the one hand, to obtain a scattering 
representation of complex systems and on the other 
hand to confirm that the scattering formalism has 
the same properties as the bond graph formalism 
and thus belonged to the formalism's class of 
network type. 

Finally, by applying the procedure described 
previously, we build the famous bond graph model 
often named: Scattering Bond Graph which will 
enable us to capture the power transfers in a simple 
and direct manner at the same time. It provides us 
with a temporal approach of the phenomena usually 
modeled with the frequencies tools. 
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Abstract 

This paper presents a detailed analysis of the 
switching loss characteristics of various existing 
discontinuous PWM (DPWM) algorithms, which use 
only one zero state and advanced DPWM (ADPWM) 
algorithms which use only one zero state with active 
state division. These DPWM algorithms reduce the 
switching losses at any operating conditions when 
compared with conventional space vector PWM 
(CSVPWM). Moreover, to reduce the complexity 
involved in the conventional space vector approach, 
the proposed PWM algorithms are developed by 
using the concept of imaginary switching times. By 
analyzing the switching loss characteristics, the 
minimum switching loss PWM algorithms are 
developed for induction motor drives. The 
theoretical evaluation is validated through the 
numerical simulation studies. 

Keywords: ADPWM, DPWM, MSLPWM, space vector 
PWM, switching loss characteristics. 

1. Nomenclature 

f s Sampling frequency 

fsw Average switching frequency 

E su b Average switching energy loss per subcycle 

V dc DC link voltage 

n a Number of switchings per sampling time 

period in phase A 
P sw Normalized switching loss 

(j) Power factor angle 

S Phase-angle 



2. Introduction 

The variable speed induction motor drives require 
variable voltage, variable frequency voltages, which 
can be obtained by using the voltage source 
inverters (VSI). Hence, nowadays the PWM 
algorithms for VSI are becoming popular. A detailed 
survey on PWM algorithms is given in [1]. Among 
the various PWM algorithms the CSVPWM 
algorithm is popular due to its numerous 
advantages compared with the sinusoidal PWM [2]. 
By utilizing the division of zero state time, several 
DPWM algorithms can be generated and a detailed 
study about this is given in [3]. Recently, novel PWM 
algorithms have been developed which use only one 
zero state same as existing DPWM algorithms, but 
involve the division of active state time in every 
sampling time interval. These ADPWM algorithms 
can be generated with space vector approach very 
easily when compared with the triangle comparison 
approach [4]. A detailed analysis of these ADPWM 
algorithms has been given in [5]. But, in the recent 
years, the research is focussed on the analysis of the 
switching loss of the inverter. Several authors have 
proposed various methods for reduced switching 
loss of the inverter [6]-[9]. However, the above 
papers deal with the conventional space vector 
approach, which involves the calculation of angle 
and sector information. Hence, the complexity 
involved in the conventional space vector approach 
is more. To overcome this problem, few PWM 
algorithms have been developed in [10] by using 
the concept of imaginary switching times. 

This paper presents differet PWM algorithm using 
the concept of imaginary switching times. The 
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Proposed algorithm reduces the switching losses of 
the inverter. Also this presets the minimum 
switching loss PWM algorithms for induction motor 
drives. 

Section 3 gives the details about the generation of 
various PWM algorithms. A detailed analysis on 
switching loss characteristics is given in section 4 for 
various PWM algorithms. Moreover, the minimum 
switching loss PWM algorithms are given in this 
section. 



3. Switching Sequences 

The proposed approach uses the instantaneous 
reference phase voltages only to calculate switching 
times. This method does not depend on the 
magnitude of the reference voltage space vector 
and its relative angle with respect to the reference 
axis. If the reference voltage vector lies in the first 
sector, then the actual switching times can be 
deduced as follows [10]: 

Ti ~V an ~^V bn =T an -T bn (1) 

v dc v dc 

And T 2 = ^f~V bn - ^f~y cn - T bn - T cn (2) 

V dc V dc 



From the above expressions, the imaginary 
switching time periods proportional to the 
instantaneous values of the reference phase 
voltages are defined as 



T = 

± an 



T = 
1 bn ~ 



v, 



v n 



dc 

\ 



v, 



V, 



T, 



dc 

\ 



bn 



( 3 ) 



\Y dc J 



V,: 



Thus, the active voltage vector switching times can 
be represented by the time difference between the 
imaginary switching time periods. The switching 
times T ari’ T bn and T cn could be negative when the 



instantaneous reference voltages are negative. 
Hence, these times are called as imaginary switching 
times. The active vector switching times can be 
calculated in each sampling interval as follows: 



TMax=Max{T an ,T bn ,T cn ) 

Let T Min =Min ( T am T bm T cn) W 

T Mid = ( T an » T bn > T cn ) 



where Max, Min and Mid are the nominal values 
used during the sampling interval. The 
function Max(T an , T bn , T cn ) selects the maximum 



value 



among 



Tan, Tbn and T cn 



Similarly 



Min(T an i T bn J cn ) selects the minimum value and 
Mid(T an ,T bn ,T cn ) selects the middle value. Finally, 



the active state times Ti and T 2 may be expressed as 

[ 10 ] 



T \ ~ T Max ' 



1 Mid 



T 2 ~ T Mid ~ T Min 



( 5 ) 



The zero voltage vectors switching time is calculated 
using (6). 



Tz = T s -T l -T 2 



( 6 ) 



The CSVPWM algorithm employs equal division of 
zero voltage vector time within a sampling time 
period. However, by utilizing the freedom of zero 
state time division, various discontinuous PWM 
(DPWM) algorithms can be generated. In the 
proposed PWM sequences the zero state time will 
be shared between two zero states as r 0 for V 0 and 



T 7 for V 7 respectively, and can be expressed as [6] 



T 0 =k o T z -,T 7 =(l-k o )T z (7) 

If k Q =0.5,0 and 1, then CSVPWM, DPWMMAX and 

DPWMMIN can be obtained respectively. 
When k Q = 0, any one of the phases is clamped to 

positive dc bus for 120 degrees over a fundamental 
interval and when k 0 = 1 , any one of the phases is 

clamped to negative dc bus for 120 degrees over a 
fundamental interval. Thus, in the first sector, 
CSVPWM uses 0127-7210 sequence, DPWMMAX 
uses 721-127 sequence and DPWMMIN uses 012- 
210 sequence. Various DPWM algorithms can be 
generated with step change of k Q between zero 

and one [3]. 

The above existing DPWM algorithms cannot switch 
more than once in every sampling time period. 
Hence, these can be generated with triangular 
comparison approach also. But in the proposed 
switching sequences, one of the phases will be 
clamped while the one of the other phases will 
switches twice in every sampling time interval. 
Hence these sequences also known as double- 
switching clamping sequences. Same as the existing 
DPWM algorithms, by changing k Q value and using 

the active sate division, advanced DPWM algorithms 
can be generated. If k Q = 0 and 1 , then 

ADPWMMAX and ADPWMMIN can be obtained 
respectively. The switching sequences pertaining to 
all six sectors for these PWM algorithms are listed in 
Table 1. Moreover, the Table 1 shows the clamping 
phase and double-switching phase also in every 
sector for ADPWMMIN and ADPWMMAX 
algorithms. The other ADPWM algorithms can be 
generated with step change of k Q between zero 

and one as given in [3]. The following Tables will 



40 



ACSE Journal, ISSN: 1687-481 1 , Volume 10, Issue 1, ICGST LLC, Delaware, USA, November 2010 



give the comparison between the existing DPWM 
and ADPWM algorithms. After examining the all 
possible advanced DPWM (ADPWM) algorithms, it 
can be observed that in each ADPWM algorithm 
though one of the phases switches twice in a 
sampling time interval the total number of 
commutations per subcycle is only three. This is 
same as the CSVPWM algorithm. Hence the 
sampling frequency of the CSVPWM and ADPWM 
algorithms is the twice the average switching 
frequency ( f s = 2 f sw ), whereas the sampling 

frequency of the existing DPWM algorithms is three 
times the average switching frequency (f s = 3 f sw ). 



4 . Switching Loss Characteristics 

The switching losses of a pulse width modulated VSI 
fed induction motor drive are load dependent and 
increase with the current magnitude. The switching 
losses of the inverter also depend on type of PWM 
method. With continuous PWM methods, all the 
three phase currents are commutated within each 
carrier cycle of a full fundamental cycle. Therefore, 
for all continuous PWM methods the switching 
losses are the same and independent of the load 
power factor angle. However, with the DPWM 
methods, the switching losses are significantly 
influenced by the type of modulation method and 
load power factor angle. 



Table 1: ADPWMMAX and ADPWMMIN sequences in all six sectors 



Sector 


ADPWMMIN 


ADPWMMAX 


sequence 


Clamping 

phase 


Double Switching 
phase 


sequence 


Clamping 

phase 


Double Switching 
phase 


I 


0121-1210 


c 


b 


7212-2127 


a 


b 


n 


0323-3230 


c 


a 


7232-2327 


b 


a 


m 


0343-3430 


a 


c 


7434-4347 


b 


c 


IV 


0545-5450 


a 


b 


7454-4547 


c 


b 


V 


0565-5650 


b 


a 


7656-6567 


c 


a 


VI 


0161-1610 


b 


c 


7616-6167 


a 


c 



Table2:Switching sequences for DPWMO and 
ADPWMO 



Sector 


DPWMO 


ADPWMO 


I 


012-210 


0121-1210 


n 


723-327 


7232-2327 


m 


034-430 


0343-3430 


IV 


745-547 


7454-4547 


V 


056-650 


0565-5650 


VI 


761-167 


7616-6167 



Table3 Switching sequences for DPWM2 and 
ADPWM2 



Sector 


DPWM2 


ADPWM2 


I 


721-127 


7212-2127 


n 


032-230 


0323-3230 


m 


743-347 


7434-4347 


IV 


054-450 


0545-5450 


V 


765-567 


7656-6567 


VI 


016-610 


0161-1610 



Table 4: Switching sequences for DPWM1 and ADPWM1: 



Sector 


DPWMl 


ADPWM 1 


0< « < 30° 


30° <a< 60° 


0 < a < 30° 


30 ° <a< 60° 


I 


721-127 


012-210 


7212-2127 


0121-1210 


n 


032-230 


723-327 


0323-3230 


7232-2327 


m 


743-347 


034-430 


7434-4347 


0343-3430 


IV 


054-450 


745-547 


0545-5450 


7454-4547 


V 


765-567 


056-650 


7656-6567 


0565-5650 


VI 


016-610 


761-167 


0161-1610 


7616-6167 
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Table 5: Switching sequences for DPWM3 and ADPWM3: 



Sector 


DPWM3 


ADPWM3 


0<a< 30° 


30 ° <a< 60° 


0 < or < 30° 


30° <a< 60° 


I 


012-210 


721-127 


0121-1210 


7212-2127 


n 


723-327 


032-230 


7232-2327 


0323-3230 


m 


034-430 


743-347 


0343-3430 


7434-4347 


IV 


745-547 


054-450 


7454-4547 


0545-5450 


V 


056-650 


765-567 
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.In the DPWM methods, the devices are clamped to 
either negative bus or positive bus for a total of 
120° and hence reduce the switching losses of 
inverter over continuous PWM methods. Hence, in 
DPWM algorithms, the load power factor and the 
modulation method together determine the time 
interval that the load current is not commutated. 
Therefore, it is necessary to derive the switching loss 
characteristics to compare the switching losses of 
various DPWM algorithms. This section presents a 
comparison of inverter switching losses due to 
conventional SVPWM and DPWM methods. 

The switching loss in an IGBT depends mainly on the 
dc link voltage (V dc ) instantaneous line current and 
turn-on and turnoff times. However, the Yd C and 
times are assumed to be constant for different 
instantaneous line currents. Hence, to study the 
switching losses of the inverter, it is sufficient to 
consider the product of instantaneous line current 
magnitude of a particular phase and the number of 
switchings per sampling time period in that phase 
( n a ), corresponding to the PWM sequence 

considered. This product is referred to as the 
switching loss factor (SLF). Since the three phases 
are symmetric, it is enough to analyze one phase 
only. The switching losses of a PWM-VSI induction 
motor drive can be modeled analytically by 
assuming linear current turn-on and turn-off 
characteristics with respect to time for the inverter 
switching devices and considering only fundamental 
component of the load current. Let the phase 
current is 

*a =/ maxSin ((Ot-</>) (8) 

where i a is the instantaneous fundamental phase 
current, 7 max is the maximum value of the 
fundamental phase current and </>is the line side 
power factor angle. The average switching energy 



loss per subcycle (E sub ) in an inverter leg is as given 



by 



„ . , 1 V a \l max sin(fflf-0)| 

E sub ( av S ) = - J 

n 0 ^max 



d(cot) 



( 9 ) 



1 71 

= — \ n a |sin (cot - (/))\ d ( cot ) 
n Q 

To obtain the measure of the inverter switching 
losses, the average switching energy loss per 
subcycle must be multiplied by the number of 
subcycles per second, i.e., the sampling frequency 
( f s ). The sampling frequency of the CSVPWM and 



ADPWM algorithms is two times the switching 
frequency ( f sw ), while it is three times the switching 

frequency (f sw ) for the existing DPWM methods. 



The average switching energy loss over a 
fundamental cycle for CSVPWM equals (2/ 7r). The 
normalized switching loss due to given PWM 
algorithm can be obtained as given in (10). 

E sub( av S )* f s 



P = 

1 sw 



fyn) Vs, 



( 10 ) 



From (10), the normalized switching loss due to 
existing DPWM methods can be obtained and given 
in (11). 

3 jr 

P sw =^E sub (avg) (11) 

By observing the DPWM0, DPWM1 and DPWM2 
modulating waveforms, it can be given that there is 
a 30° phase-angle ( 5 ) distance between their de- 
link clamped 60° segments. Hence, a new PWM 
algorithm can be introduced as generalized DPWM 
(GDPWM) algorithm which covers the DPWM0, 
DPWM1 and DPWM2 algorithms. The 8 and (j) 



dependent switching phase current and normalized 
switching energy loss per subcycle waveforms of 
GDPWM algorithm is shown in Fig.l. The GDPWM 
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algorithm can be represented by using space vector 
as shown in Figure 2. For 8 =0°, 30° and 60°, 
DPWM0, DPWM1 and DPWM2 algorithms can be 
obtained respectively. The variation of average 
switching energy loss per subcycle ( E sub ) over a 

fundamental cycle for various DPWM algorithms for 
different values of 8 and $ is shown in Figure 3- 

Figure 6. From these, it can be observed that at 
unity power factor DPWM1 clamps a phase around 
its current peak, which leads to a significant 
reduction in E sub (avg) over the remaining DPWM 

and conventional SVPWM algorithms. Therefore, the 
switching loss mainly depends on the power factor 
angle by which the line current lags/leads the line 
voltage. Moreover, it may be observed that the 
effect of different PWM sequences on the switching 
loss does not depend on V re f or the fundamental 

frequency. 





Figure 1: the average switching loss of GDPWM 




Figure 2: space 
algorithm 



vector illustration 



of GDPWM 





Figure 3(a) 

From Figure 3 - Figure 6, the expressions for 
normalized inverter switching loss corresponding to 
different DPWM algorithms can be derived as given 
in [3]. The variation of normalized switching loss of 
different DPWM algorithms with power factor angle 
(both at lagging and leading) is shown in Figure 7. 





Figure 3(b) 

Figure 3. Variation of normalized switching energy loss for 
DPWM0 method: (a) (f) = 0° (b) </> = 30° 





Figure 4(a) 
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Figure 4(b) 

Figure 4: Variation of normalized switching energy loss for 
DPWMl method: (a) </> = 0° (b) </> = 30° 





Figure 5 (a) 





Figure 5 (b) 

Figure 5: Variation of normalized switching energy loss for 
DPWM2 method: (a) (/) = 0° (b) (/) = 30° 
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Figure 6 (a) 





Figure 6 (b) 

Figure 6: Variation of normalized switching energy loss for 
DPWM3 method: (a) <j> = 0° (b) </> = 30° 

From Fig. 7, it can be found that the DPWMl gives 
minimum switching loss at unity power factor, 
because DPWMl clamps at peak phase current at 
unity power factor. Similarly, DPWM0 and DPWM2 
give minimum switching losses at 30° leading and 
lagging power factors respectively. But, DPWM3 
gives minimum switching loss near the zero power 
factor. The value of 8 can be varied in accordance 
with (j) to achieve reduction in switching losses. 




Figure 7 The variation of normalized switching loss of different 
DPWM algorithms (A: DPWM0; B: DPWMl; C: DPWM2 and D: 
DPWM3) 



The minimum switching power loss solution of the 
DPWM0, DPWMl and DPWM2 yields to GDPWM 
algorithm. The optimal solution of GDPWM is 
obtained by selecting 8 = (;r/6)+ $ 

for-(;r/6)< £ <(;r/6). Outside this range 8 = 60° for 

lagging power factor angle and 8 = 0° for leading 
power factor angle. The normalized switching loss 
variations for different DPWM algorithms and 
GDPWM algorithm are shown in Fig. 8. This shows 
that the GDPWM algorithm combines the DPWM0 
and DPWMl and DPWM2 algorithms and also gives 
less switching loss.The GDPWM does not consist of 
DPWM3 algorithm. But from Fig. 7, it can be 
observed that the DPWM3 gives minimum switching 
losses outside the — (5 tt/12) < $ < (5tt/ 12) range. In 
this paper, by considering DPWM3, the minimum 
switching loss PWM (MSLPWM) algorithm with 
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existing DPWM algorithms is proposed. The 
variations of normalized switching losses of various 
DPWM algorithms and proposed minimum 
switching loss PWM (MSLPWM) algorithm are 
shown in Figure 9. The comparison between 
GDPWM and proposed MSLPWM is shown in 
Figure 10. 




Figure 8: The variation of normalized switching loss of different 
DPWM algorithms with power factor angle (A: DPWMO; B: 
DPWM1; C: DPWM2; D: DPWM3 and E=GDPWM) 



Similar to the existing DPWM algorithms, the 
switching loss analysis of ADPWM algorithms can be 
carried out. The normalized switching loss due to 
ADPWM algorithm can be obtained same as 
previous chapter and given by 

TC 

P sw, ADPWM = — E sub( av 8 ) (12) 




Figure 9: The variation of normalized switching loss of different 
DPWM algorithms with power factor angle (A: DPWMO; B: 
DPWM1; C: DPWM2; D: DPWM3 and E= MSLPWM) 




Power factor angle (degrees) 

Figure 10: The variation of normalized switching loss of GDPWM 
(A) and MSLPWM (B) algorithms 



The modulating waveforms of ADPWM0, ADPWM 1, 
ADPWM2 and ADPWM3 are similar to 
corresponding existing DPWM algorithms. Hence, 
same as existing DPWMO, DPWM1 and DPWM2 
algorithms, there is a 30° phase-angle (^) distance 
between their dc-link clamped 60° segments in the 
ADPWM0, ADPWM1 and ADPWM2 modulating 
waveforms also. Hence, a new PWM algorithm can 
be introduced as advanced generalized DPWM 
(AGDPWM) algorithm which covers the ADPWM0, 
ADPWM 1 and ADPWM2 algorithms. The variation of 
average switching energy loss per subcycle ( E sub ) 

over a fundamental cycle for various ADPWM 
algorithms for different values of 8 and $ is shown 

in Figure 11 - Figure 14. 





Figure 11 (a) 




Figure 1 1 (b) 

Figure 11: Variation of normalized switching energy loss for 
ADPWM0 method: (a) (/) = 0° (b) (/) = 30° 



At the same time ADPWM3 is not effective in 
reducing the switching energy loss at unity power 
factor. At zero power factor lagging, ADPWM3 
reduces the E sub (avg ) . Thus, both in existing as well 



as advanced DPWM algorithms, DPWM1 and 
ADPWM1 are better in terms of switching losses at 
power factors close to unity, while DPWM3 and 
ADPWM3 are better at power factors close to zero. 
From Figure 11 - Figurel4, the expressions for 
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normalized inverter switching loss corresponding to 
different ADPWM algorithms can be derived as 
given in [8]. The variation of normalized switching 
loss of different ADPWM algorithms with power 
factor angle (both at lagging and leading) is shown 
in Figure 15. 




Angle(degrees) 



Figure 12 (a) 




Figure 12 (b) 

Figure 12: Variation of normalized switching energy loss for 
ADPWM 1 method: (a) (j) = 0 ° (b) <j) = 30 ° 






Figure 13 (b) 

Figure 13. Variation of normalized switching energy loss for 
ADPWM2 method: (a) (/) = 0 ° (b) </) = 30^ 





Figure 14 (a) 

Figure 14. Variation of normalized switching energy loss for 
ADPWM3 method: (a) (/) = 0 ° (b) (/> = 30° 




Figure 15: The variation of normalized switching loss of different 
ADPWM algorithms with power factor angle (A: ADPWM0; B: 
ADPWM 1; C: ADPWM2 and D: ADPWM3) 
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From Figure 15, it can be observed that the 
ADPWM1 gives minimum switching loss at unity 
power factor, because ADPWM1 clamps at peak 
phase current at unity power factor. Similarly, 
ADPWMO and ADPWM2 give minimum switching 
losses at 30° leading and lagging power factors 
respectively. But, ADPWM3 gives minimum 
switching loss near the zero power factors. The 
value of 8 can be varied in accordance with $ to 

achieve reduction in switching losses. The minimum 
switching power loss solution of the ADPWMO, 
ADPWM1 and ADPWM2 yields to AGDPWM 
algorithm. The optimal solution of AGDPWM is 
obtained by selecting 8 = { 71 / 6)+^ for 

-(tt/6)< 8 < {tt/6) . Outside this range 8 = 60° for 

lagging power factor angle and 8 = 0° for leading 
power factor angle. The normalized switching loss 
variations for different ADPWM algorithms and 
AGDPWM algorithm are shown in Fig. 16. This 
shows that the AGDPWM algorithm combines the 
ADPWMO and ADPWM1 and ADPWM2 algorithms 
and also gives less switching loss. 




Figure 16: Variation of normalized switching loss of different 
ADPWM algorithms with power factor angle (A: ADPWMO; B: 
ADPWM 1 ; C: ADPWM2; D: ADPWM 3 and E=AGDPWM) 



The AGDPWM does not consist of ADPWM3 
algorithm. But from Fig. 15, it can be observed that 
the ADPWM3 gives minimum switching losses 
outside the — (5 tt/12) < $ < (5;r/12) range. In this 

papaer, by considering ADPWM3 also into account 
advanced minimum switching loss PWM 
(AMSLPWM) algorithm with ADPWM algorithms is 
proposed. The variations of normalized switching 
losses of various ADPWM algorithms and proposed 
AMSLPWM algorithm are shown in Fig. 17. The 
comparison between AMSLPWM and MSLPWM is 
shown in Fig. 18. From Fig 18, it can be observed 
that the proposed AMSLPWM algorithm gives 
minimum switching loss over the MSLPWM in the 
- G/6) < (Z> <(;r/6) range. At the same time, 
MSLPWM algorithm gives minimum switching 
losses outside the - (71/6) <</> <(71/6) range. 




Figure 17: variation of normalized switching loss of different 
ADPWM algorithms with power factor angle (A: ADPWMO; B: 
ADPWM 1 ; C: ADPWM2; D: ADPWM 3 and E= AMSLPWM) 




Figure 18: The variation of normalized switching loss of MSLPWM 
(A) and AMSLPWM (B) algorithms 

5. Conclusions 

Space vector based minimum switching loss PWM 
algorithms are presented in this paper. The proposed 
algorithms did not use the information regarding sector 
and angle calculation of switching times and hence 
reduces the complexity involved. The minimum 
switching loss PWM algorithms are derived using 
existing DPWM and ADPWM algorithms and compared. 
From the results, it can be observed that AMSLPWM 
algorithm gives minimum switching loss over the 
MSLPWM in the - {ji/ 6) <<j>< (tt/6) range. At the 
same time, MSLPWM algorithm gives minimum 
switching losses outside the - (;r/6) < (j) < (71/6) 

range. 
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Abstract 

In this paper, we present an adaptive fuzzy control 
for a class of unknown nonlinear systems. The 
strategy of control uses the adaptive fuzzy system 
Takagi-Sugeno (T-S) to approximate the part of the 
primary control. In order to guarantee the stability 
and high performance, the auxiliary part is 
incorporated in the control law. The proposed 
compensation control is combined a Hermite 
function and sliding mode control (SMC). The 
Hermite function is used to predict and to remove 
the fuzzy approximation and attenuates the 
influence of the external disturbances. However, the 
SMC is designed to ensure the stability and a good 
tracking. The parameters of the control are adjusted 
on line by the adaptive law with stability and 
convergence analysis using the Lyapunov approach. 
The robust and the stability of the control scheme 
are proved and simulation results are given to verify 
the effectiveness of the proposed approach. 

Keywords: Adaptive control, Hemite function, Sliding 
mode, Fuzzy logic, Nonlinear system 

1. Introduction 

The control of nonlinear systems has been an 
important research topic [1, 2]. Traditionally, control 
system design has been tackled using mathematical 
models derived from physical laws. In fact, most of 
the parameters and structure of the system are 
unknown due to environment changes, modelling 
error and unmodelled dynamics. Then, the designer 
is confronted with the above problems when 



designing control systems. Face the above 
problems; there exist several techniques, in 
particular, the intelligent technologies as neural 
networks, fuzzy logic, genetics algorithms and 
evolutive computation [3, 9, 10]. 

Fuzzy logic, as one of the most useful approaches 
for utilizing expert knowledge, has been an active 
field of research during the past decade [9, 11]. 
Fuzzy logic control has found promising 

applications for a wide variety of industrial systems 
specifically applicable to plants that are 

mathematically poorly modelled [7, 16]. 

The ability of converting linguistic descriptions into 
automatic control strategy makes it a practical and 
promising alternative to the classical control scheme 
for achieving control of complex nonlinear systems. 
Recently, a great amount of the effort has been 
devoted to describing nonlinear system using a T-S 
fuzzy model. Based on the universal approximation 
capability, many effective adaptive fuzzy control 
schemes have been developed to incorporate with 
human expert knowledge information in a 
systematic way, which can also guarantee stability 
and performance criteria [3, 4, 10]. 

There have been many successful applications in 
fuzzy control in recent years. In spite of the success, 
there are still many basic issues that remain to be 
further addressed. Stability analysis and systematic 
design are certainly among the most important 
issues for fuzzy control systems [14]. In this respect, 
several direct and indirect adaptive fuzzy control 
schemes have been introduced for controlling 
nonlinear systems [15]. In the direct scheme, the 
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fuzzy system is used to approximate an unknown 
ideal controller. On the other hand, the indirect 
scheme uses fuzzy systems to estimate the plant 
dynamics and then synthesizes a control law based 
on these estimates. In these adaptive fuzzy control 
schemes, the controllers are generally composed of 
two main components. 

This paper proposes an adaptive fuzzy control for 
nonlinear system with unknown dynamics. The 
strategy of control is based on fuzzy system 
combined Hermite function and SMC. The primary 
control action is approximated by the adaptive fuzzy 
model type T-S. The auxiliary part of control is 
incorporated to ensure the stability and the tracking 
performance. The design of this compensation 
control is based on Hermite function and SMC. The 
fuzzy approximation error and the influence of the 
external disturbances are predicted by using the 
orthogonal Hermite basis function. The SMC is 
designed to guarantee the stability and robustness 
of the proposed approach. The fuzzy and Hermite 
basis parameters are estimated on line by the 
adaptive laws with stability and convergence are 
analysed in the Lyapunov sense. 

The paper is organized as follows: the problem 
formulation is presented in section 2. In section 3, 
the adaptive fuzzy control is proposed. Simulation 
results are presented in section 4. 



2. Problem Formulation 

Consider a nonlinear system described as follows: 
x = 



x 2 = x 3 



x n = f(x 9 t) + bu + d(t) 
y = x i 

Or equivalent: 

J x (n) = f(x,t) + bu + d(t) 
where x - [v. x, V" 



( 1 ) 



( 2 ) 



2 ’ • 



■>* n 



f sR' 



the state vector of the systems which is assumed 
measurable. The function / is unknown and 
nonlinear, and b is constant positive, u e R and 
y e R are respectively the input and output 
system. d(t ) is the unknown external disturbance. 
Our objective is to design an adaptive control which 
will cause the output y to track a desired reference 

y r (i.e., e 0 — > 0 ) and have up to (n - 1 ) bounded 
derivatives. 

The tracking error defined by: 



e 0 = y - y r 0 ) 

The control law is designed to have the following 
idealized control [3]: 

u =7 ~(-f(x,t) + yl n) (4) 
b 

It is desired that the tracking error of the system 
follow: 

e ( 0 n) + k n _ x e { ;~ X) + • • • + k 0 e 0 = 0 , where 

k { (i = 1 , 2 ,...., ft — 1 ) are the coefficients of the 
Hurwitz polynomial: 

H(z) = z n ~ l + k n _ x z n ~ 2 +--- + k 1 . 

In fact, most the parameters and structure of the 
systems are unknown due to environment changes, 
modelling errors and unmodelled dynamics. Thus, 
the primary control law (4) can not be implemented. 
To solve this problem, the direct adaptive fuzzy 
control scheme is proposed to approximate the 
idealized control (4). 

3. Direct adaptive Fuzzy Control 

In this section, direct adaptive fuzzy control is 
designed, with guaranteed stability of integrated 
closed loop system. The direct adaptive controller is 
designed as follows: 

u = u(x,6) + u h + u s (5) 

where u(x,0) is the direct adaptive fuzzy control, 
the auxiliary control u h is designed to predict the 
fuzzy approximation error and the influence of the 
external disturbances, the control u s is incorporated 

in the control law in order to ensure the stability 
and a good tracking error. 

The adaptive fuzzy system T-S is used to 
approximate the primary control action u . In the 
respect, the parameters of the fuzzy control is 
updated on-line using the Lyapunov approach to 
guarantee the stability and convergence analysis. 

The T-S fuzzy model is described by fuzzy if-then 
rules, which locally represent linear input-output 
relations of a system. The fuzzy system is of the 
following form: 

R l :If x x is F[ and x 2 is F[ . . m ,x n is F l n 

( 6 ) 

then u l is 0 l 

l = 1,2 

where R l represents the Ith fuzzy inference rule, 
x i ,i = l are measurable variables of system, 

F* denotes fuzzy set. Each fuzzy set F* is 
associated with a membership function ju f i and the 
center of this membership function is the operating 
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point. By using the weight average fuzzy inference 
approach [12], we obtain: 



u(x, 0) 



Zfl »F} 0 1 



1=1 i = 1 



zft^ 

1=1 i=l 



(7) 



= 0 T Z(x) 

0 = [0 l ,0 2 ,...,0 m ] T is parameter vector and 

£(*) = [^(x),^ 2 (x),...,^ m (x)] T \s a regressive 

vector with regressor £,(x)(( 1 <1 <m ) , (m is the 

number of rules), which is defined as fuzzy basis 
function: 



FK 

ith 

Z=1 i=l 

The adjustable fuzzy parameters of w(x,#)are 



tuned on line using the Lyapunov approach. In order 
to guarantee that the parameters are bounded, we 
introduce the projection algorithm [12] to restrict 
them in the closed set Q . 



- y x ab^(x) if (|/9| < M or 



d — < 



(|#| =Mand y l abO T ^(x)> 0)) 



Projf-^ crb^(x)] otherwise 



(9) 



a(x,t) = 0 represents a time varying sliding 

surface which is defined in the state space/?” by the 
scalar equation: 

cr(x,t) = e q” 1} + k n _ x e q” 2) 4 t- (10) 

where ^is fixed adaptive gain. Proj[J represents the 
projection operator which is defined by: 

Proj \-y l a b<^(x)\ = -y x ab<!;(x) + 



functions Hermite. This method consists on 
projecting w = u-u*ar\d d. 

The well known continuous Hermite functions are a 
basis for the set of polynomials, defined on time 
domain ]- oo, + oo[, and are given by: 

(-l)*exp(F) 

S k (t) = -f —yexp(-t 2 ) (12) 

VFFVT dt 

Where k = 0,1,2... is the function order. These 
functions can be written in the vector form by 
defining q the order state vector: 

S(t) = [s 0 (t) 5,(0 ... S q (t)\ 



The value at t of any Hermite polynomial S k can be 

determined using the following 

recursion (0 <k<q)\ 

Where: 

1 /o^- -f ^ 

5 0 (0 = ^exp 5,(0 = 7^exp(--) 

y 7i i yn i 

Then the control law u h is derived and based on 

the Hermite function basis to attenuate the 
approximation and the influence of the 

disturbances. 

u h =-(pS T {t ) (13) 



The adjustable parameter vector: 

<p = [cp x , (p 2 ,...,(p q ] T is estimated on line using the 



Lyapunov approach: 

-y 2 crb 5(0 if (|<^ <(Vor 



cp=< 






=Nandy 2 crb Sit )>())) 



Pr°j[-r 2 abS(t)] otherwif 



(14) 



y.crb 



6d T %(x) 



(11) 



In the practice, however, usually not all states are 
available, and the hypothesis which consists on 
considering that the nonlinear dynamics to have the 
behaviour of the linear combination of known 
nonlinear functions seems more restrictive, since 
generally the nonlinear dynamics are completely 
unknown. Therefore, the prediction of the later is 
which a great amount of interest in the control 
nonlinear adaptive control. The fuzzy approximation 
error and the influence of the external disturbances 
are estimated by using the orthogonal basis 



The strategy of control u s is based on SMC 

approach to ensure stability, tracking and consistent 
performance: 

u s = -k d sat(a(xj)) (15) 



Theorem: Consider the nonlinear system (2). If the 
direct adaptive fuzzy control (5) is applied with the 
parameters vector 6 and cp are adjusted by the 
adaptive law (9)-(14). The closed-loop system is 
stable in sense that all the signals are bounded and 
the tracking performance is guaranteed. 



Proof: 

Consider the following Lyapunov function as: 
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V = -a 2 + — <M> r + — 

2 2 y x ly 2 

where 

O = 0* - 0 and = cp * —q> 

The optimal parameters are given by: 
6* = argmin(sup x&R « u* -m(x,^)|) 



(16) 






(p = arg rmn(sup re/? » 

<P^h 

where Q u = {#/|#| < m} and 
Cl h = \<pl\(p\ < Afjare the convex compact sets 
which contain feasible parameter sets for# and cp . 
M and N are pre-specified constants, 
w =b(u(x 9 0*) - u*) + d (19) 

The time derivative of V (16) is: 

V = cr& + — OO r + — 



w- 





(17) 


■PS T ) 


(18) 



( 20 ) 



r i r 2 

*\n - IL 



& = x (n) - y ( r n) + 7;_,M (|) = b(u(x,0)-u*) 



( 21 ) 



+ Z?W ;t +d(t) 

V = <7 b(u(x,0) - u* ) - <7 bu h + crbu s + 

1 . I . (22) 

crj(0 + — OO r + — v F v F r 

Yx Yi 

V = crb(u(x,0)-u(x,8*) + u(x,0*)-u*) 



-crbu h +crbu s +crd(t) + —(bO T + — 't*'? 7 

Yx Yi 

V = crb(w + d -u h +u s ) + crb(u(x,0*) 

1 . 1 . (23) 

-u*) + — OO r + — 'F'F 7 

Yx Yi 

V =crb(w + d-u* h +u* h -u h +u s ) 

+ ab(u(x,0*)-u*) + + — (b® T + — x ¥'¥ T 

Yx Yi 

Define the minimum approximation error: 

£ = w + d -u* h (24) 

The V can be written as: 

V = <7 b{s + u s ) + <7 b(u* h ~u h ) + 

(7b{u{x,0*)-u ) + + — <& <D r + — v F v F r 
Yx Yi 

V = ab(£ + u s ) + ab x ¥*S T +crM>*£ r + 



— OO r + — 



(25) 



Yx 



Y 2 



By consideration of the update laws (9)-(14): 

V = crb(£ + u s ) + — (6 + (7b%(x)) + 

Yx 

— ^ T (^ + c7bS) 

Y 2 

Hence 

V <\<7\(\e\-k d <7) 

So 

V <0 



( 26 ) 



(27) 



4. Simulation Results 

In this section, we test the proposed adaptive fuzzy 
control on the regulation control of the 
servomechanism. The dynamic equations of the 
servomechanism can be described in space state as 
[15]: 

JCj = x 2 

< x 2 = f(x,t) + u + d(t) (28) 

7 = *i 

/ (x, t) = -x 2 - 0.4 sin(x 1 ) . 

In the first step, we need to define some fuzzy sets 
to cover the state space. The choice of the number 
of fuzzy set and the constant M is related to 
knowledge of expert on the system. We consider 
M =16 and three fuzzy membership functions are 
chosen as in figure 1: 




Then there are 9 rules to approximate the primary 
control law/2(x, #) . 

The fuzzy rules are defined by the following 
linguistics description: 

R l :If x 1 is F[ and x 2 is F 2 then u l is 0 l (29) 
Where / =1,2,3 and m = 9 is the number of the 
fuzzy rules. 
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The choice of the number of fuzzy set and the 
constant M is related to knowledge of expert on the 
system. The control objective is to maintain the 
system to track the desired angle 

trajectory: * (8111(0.0250 . 

The sliding surface is defined as: 
cr(x, t) = k x e 0 (t) + k 2 e 0 (t) 

with k x =5 and k 2 =l. The initial consequent 
parameters 0(0) of fuzzy control and ^(0)are 
respectively chosen randomly in the interval [-1, 1] 
and [-3, 3]. The fuzzy approximation error and the 
influence of the external disturbances are predicted 
using the orthogonal basis of Hermite function 

S(t) = [s 0 (t) 5,(0 5 2 ( 0 ]. 




Figure 2: Reponses of y(t) and y r if) 




Figure 3: The control signal 

Let the learning rat ey x = 0.5, y 2 = 0.3 . From figure 
2, it can be seen that, the tracking performance is 
obtained with unknown nonlinear dynamics and in 
presence of disturbances. The corresponding fuzzy 



control signal is given in figure 3. The figures 5 and 
6, illustrates the behaviour and the convergence of 
the parameters of fuzzy control and Hermite 
function. 
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Figure 4: Evolution of 0 



0 20 40 60 30 100 120 140 160 180 200 
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Figure 5: Evolution of Cp 



5. Conclusion 



In this paper, an adaptive fuzzy control has been 
proposed for a class of unknown nonlinear systems. 
We introduced the fuzzy system T-S to approximate 
the part of the primary control. The auxiliary part of 
control is composed on two terms. The first is 
incorporated to predict the fuzzy approximation 
error and the influence of the external disturbances. 
The synthesis of the second is based on the SMC 
approach in order to ensure the stability and the 
tracking performance. The closed loop system is 
stable in the sense of Lyapunov approach. The 
simulation results illustrate the effectiveness and 
robustness of the proposed fuzzy control method. 
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